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A. BACKGROUND 


For a few years the Surface Effect Ship (SES) has been 
receiving increasing attention in the United States Navy 
(USN), and detailed studies have been undertaken since then. 
The two major categories of SES are tne Air Cushion Venicles 
(ACV) and the Captureä Air Bubble Ships (CAB). 


Whereas the weight of the ACV is entirely supported гу 
the pressure differential in the plenum chamber all 
surrounded by flexibie seals between the water surface and 
the hull structure, the CAB uses rigid sidewalls extending 
Шо спе water and thus giving additional support to that of 
the captured air bubbie. Some lift forces are nade up by the 
flexible pow and stern seais as well as aerodynamic force 


at higher speeds are usually taken into account. 


The latter type of craft has become of primary interest 
to the USN for ocean going applications and led to the 
construction of several test crafts, among those the 
Aerojet-General 100-A, the Bell Aerospace Systems 100-8, 
both of about 100 tons displacement, and the XR-3 with 
approximately 3 tons displacement. This study will be 


concerned with the XR-3. 


This ship was built by the Naval Ship Research ana 
Development Center (NSRDC) in 1966 and tested for several 
years. In March 1970, the XR-3 was delivered to the Naval 
e graduate School (NPS), Monterey, California, for further 
investigations. After a digital conputer simulation program 
for the 100-B test craft was developed by Oceanics,Inc., leo 
and Boncal (Ref. Т] converted this program to represent the 
XR-3 test craft as there were substantial design differences 


between those two ships. 
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In December 1974, Forbes (Ref. 2] undertook a validation 
of the Loads and Motion computer program in calm water and 
introduced some nodifications үп various program 
Subroutines. At the same time Finley [Ref. 3] refined the 
bow and stern seal subroutines and developed new fan maps in 
order to achieve more realistic air flow rates and thus an 


inproved craft performance. 


DODJECTIVES 


The purpcse of this thesis is the study of the pitch and 
roll response of the XR-3 Loads and Motions computer program 
in a comparison between the two modified versions of Refs. 2 
and 3. Changes are introduced, where neccessary, with the 
idea to improve the computed pitch and roll behaviour. The 
final goal is to select an optinal model from the results 


obtarned for future studies. 


In order to ease the further discussion the computer 


models involved are named as follows: 


ama 


Progran 1 computer simulation program used in Ref. 


t 


Program 2 computer simulation program used in Ref. 


Program 3 computer simulation program used in Ref. 3 

A users manual for the final program chosen is attached 
in Appendix A and shoulä be read in conjunction with this 
text in order to better understand the subroutine programs 


and variables under discussion. 
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EI. GENET ALT DISCUSSION 


ee чьашь сне чыдыг ш =: = un ma — اف سک‎  — سی‎ ee сч 


A. ANALYSIS CF MOMENTS 


The development of the equations of motion considerino 
the craft to have six degrees of freedom wiil not be 
discussed here as it is well presented in Refs. 1 and 4. 
However, the forces and moments involved need some 
consideration. Іп general , various elements contribute 
separate forces to the different degrees of freedom and can 


be represented by: 


Е - Е + F + F + F 
bowl seal stern seal sidewalls rudder 
+ F + F + F +F 
propulsion waves aerodynamic airbubble 
4- 
gravitational 
With respect to the local reference system the x- and 


z-components of the individual forces determine the moments 
about the lateral y-axis and, therefore, the final pitch 
angle at steady state. These moments are calculated in the 
various subroutine programs (see Appendix A) and then added 
up to give the net moment which eventually reduced to zero 


in calm water steady-state. 


In figure 1 these moments for the XR-3 are shown in a 
schematic manner to indicate their predominant effect on the 
pitch behaviour. Counterclockwise acting moments are 


considered positive by convention. 


Іп а sirilar fashion, the moments around the x-axis 
determine the roll angle, which - provided the craft is 
synmetric about its longitudinal axis including weicht 
distributions - tends to be zero in calm water. In order to 
give an impression of the magnitude ratios the rounded 
moment values for a typical straight run at 20 knots are 


a 


given in table í. These figures were obtained from program 2 
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ENS original versionun.(See also D.1.a.) 


Table I 


= 


سے ニニ‏ ے ہے ^o AA AU Q m s we‏ سے سے = حت ہے جس یی 


(Program 2, original version) 


Component Positive Negative 

Name Moment Honent 

(lb-ft) (lb-ft) 
FHBS 705 Be 
FMSS ーーー 112 
FMSW ーーー 3827 
FARUD тт 85 
FMP 1147 ーーー 
ГЕДЕР 645 ЕСЕ 
FMBUB 1805 --- 
FWAVZ == ۳ئ‎ 
1 ۷ 0 дана ке 
Summary 4302 4302 


B. INITIAL CONDITIONS 


Reference 2 points out that the choice of the program 
ШШІЕТаі conditions is somewhat critical. If they are not 
close to the steady state values the imbalance in forces 
generated within the program will cause an execution stop 
due to an improper integration step size. Thus a tedious 
method with elaborate curves families was used in Ref. 2 to 
develop the correct data including a proper integration step 


өлде, 


References 2 and 3, unfortunately, give no information 
about the effect of their program changes on the initial 
conditions nor are the values used for the simulation runs 


completely listed. Refference 2, however, presents sone 
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steady state values for straight runs with pitch angles 
considerably larger than those reported by Olmstead[Ref. 5] 
closer to the values observed with the present 


which come 


craft configuration. 


As further data were not readily available, the initial 
conditions given in Ref. 1 were chosen and are 
table II. 


in both cases that the 


repeated in 
The first trial runs with programs 2 and 3 showed 


disturbances Mtroduced were well 


within the limits of the programs' capability. Steady state 
was reached after an avarage of less than ten seconds 
Simulated run time of the craft. 
Table II 
Initial Conditions 
Speed Pitch Draft Pienun Thrust Rudder 
angle pressure 
(knots) (deg) (in) (psf) (ib) (деа) 
10.0 1: 72 6.60 24.86 504.0 0.0 
1225 1.31 6.50 24.86 420.0 : 
19.0 0594 6.35 24.86 39220 0.0 
[5 0.60 0.15 24.86 Bog 
20.0 0:29 6 2. 86 424.0 20 
2249 0.13 5. 78 24.86 464.0 
25.0 0 09 5. 50 раб 54250 
21.5 0.06 5.30 24.86 568,0 0. 
30.0 0.05 3.02 24.86 626.0 : 


EXEC RAFT ГАТА 


The following data were used throughout the entire 


unless 


into the program (subroutine INCON). These names are used in 


this study without further explanation which can 


Otherwise 


Specified. 


The 


пишаа сасе the block numbers under which the values 


numbers in 


parentheses 


are 


be 





| 


up in Appendix A. 


XS (2) = 10.05 feet forward of transon 
SOT (3) = 24.7 feet 

XLSW (4) - 21.9 feet 

XSSI(5) = 3.97 Feet forward of transon 
XLF (5) E 3562667 

X BST (6) = 23.u4 feet forward of transon 
XBF (6) = 4.06 feet 

XL5W (7) = 200 22: 

ХРЕҮ (7) - 33 2 feet forfard of transom 
XL (7) = 20.0 feet 

XCPO (7) = 10.4 feet forward of transom 
X PO (8) =" 1.275 feet forvard of transom 
XRO (9) = -1.125 feet forward of transom 
XLAERO(10)- 20.0 feet 

СЕВАН (9) = 1.٦215۶6۰۳ 

YSW (4) = 5.37 feet fron centerline 
AVEMSW (4) = 0.5 feet 

WIDTH(7) = 10.0 feet 

YPO (8) = 5.55 feet trom centerline 
YR(9) - 5.55 feet from centerline 
BEAM(10) = 10.0 feet 

it 2) = 2.54 feet above keel 

551 (5) = 1.875 feet above keel 
ELMAXS(5) - 1.875 feet (Program 3 only) 
ZBSI(6) = 1.875 feet above keel 


ELMAXB (6) 
СЕНСАВ (6) 


1.875 feet (Program 3 only) 
1.875 feet (Program 3 only) 


! 


XBBW (7) = 10.0 feet 
BUBHGT (7) = 1.915 feet 
ZPO (8) = -0.604 feet above keel 
Z RO (9) = 70.208 feet above keel 


WEIGHT (2) =6050.0 pounds 
ШӘ 5) = 6 344 degrees (Progran 2 only) 


26 





TEBSI(6) = 63.4 degrees (Program 2 only) 





АХУ (2) -2870.0 slug-square feet 

AIYY (2) =9320.0 Slug-square feet 

А122 (2) =10580.0 slug-square feet 
ATXZ(2) =-2800.0 slug-square feet 

ШИВАК (5) = 26.0 square feet 

DPSS (5) = 1.0 pounds per square foot 
DPBS (6) = 1.9 pounds per square foot 
BLEAK (6) = 0.1 Square feet 

RAREA(9) = 0.68 square feet 

CFSW (4) г O. 

CDS (ц) = 1.28 

e255(5) = 0.9 

CFBS(6) = 0.9 

FICHRITI(7, = 0.556 

RASPR(9) = 25 

RTC (9) SI uS 

STHS (16) = 0.01 

STHP (16) = oM 


2] 


III. PITCH ANGLE 
(Steady State Conditions) 


A. PROGRAM 2 


1. Initial Runs and Evaluation 


ーー шшк у с — — ei سب — حسمب سی‎ wa c s aoe سوج بح‎ 


Using the initial conditions from Ref. 1 the steady 
State values of program 2 were obtained as reported in table 
атас 


lay in the expected range, the pitch angle increased 


III. Whereas the values for plenum pressure and 
thrust 
with higher speed. This tendency was the first point for a 


Dore thorcugh investigation of the changes made in Ref. 2. 


En VA Q fe > E A A a r s m O q A س‎ se = 


(Program 2, original version) 


Speed Pitch Draft Plenum Thrust Rudder 
angie pressure 

mots) (deg) (in) (psi) (1b) (deg) 
10.0 1.74 8.18 ZI 405.16 
025 1.71 8.10 24.46 356.42 0.0 
ПОО 1:92 7.97 24.46 3 ۶ 
17.5 1.80 7.78 24.46 373.96 
20.0 1.90 7.54 24.46 agen esos) 
Оос 201 1223 24.46 457.84 0.0 
ро 0 2.15 6.87 24.46 DON 7٣٦ 

2775 2.29 6.45 24.46 560.98 
30.0 2.47 6.01 24.33 613.63 0.0 


The simulation derived in Ref. 


contained a 


speed 
dependent correction term for the x-coordinate of tne center 
with 


of pressure which was assumed to move aft increasing 


e aana thus reducing the pitch angle [Ref. 6]. This tern 


had been removed in Ref. 2 with the assumption that due to 


28 








the existance of a pressure wave (Ref. 7] an incremental 
wedge of volume has to be added to the plenum volume which 
Varies with speed and moves the center of volume aft with 
increasing speed. The according changes introduced in 
Subroutines RHS and NAVES were expected to take care of the 


eorreet pitch angles. 


Jn order to study the effect of the previously 
mentioned change - here called water slope correction - the 
moments of the original program 2 were recorded as shown in 
table I. Then, one after the other, the water slope 
@eemeection terms in the subroutines INCGN, BOWSL, STNSL, &HS 
and SIDEWL were removed. Due to tne nature of these terms a 
different pitch angle couid be expected only at lower speed. 
In all cases, at 10 knots for example, tne angle was lower 
than the original value. The biggest change in the 
corresponding moment could be observed in the bow seal 


routine when the pitch angle dropped to 1.21 deg. 


The overall results iron tnese first trial runs 
indicated that the water slope correction provided a 
neccessary curvature in the pitch curve over the speed range 
of interest pointed out in the following investigation, but 
could not take care of declining angles with increasing 


Speed. 


2. Selection of Pitch àngle Curve 
AS pointed out in Rer. 5 the actual pitch angles of 
the XR-3 are dependent on the loading condition of the craft 
and the seal adjustments. Several curves are presented from 
which the fcllowing two were chosen: 
-—The upper limit for the pitch angles was Set by the 
so-called "ideal" curve with theoretically derived values. 
- The lower linit- was given by the so-called “three turns" 


curve. These values were actually measured on the craft 
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after the seals were adjusted with "three turns" by 


of a spindle. 


The values for these two curves 


angles of program 2 are shown 


means 


and the original pitch 


in table IV and grapheä in 


figure 2. Additionally the values needed later as initial 
conditions were plotted using subroutine COLFIL after 
Ming for thrust as a variable (plots 1 through Ч). 


Fron 


discussion was obvious. 


то m = ہے‎ 


these data the neccessity to correct the program under 


> سے س س سے سے ج کے چ‎ e em 


(Original version) 


Proqram 2 


Speed X l"Ideai" X "Three Turns" Program 2 
(knots) (deg) (deg) (deg) 
1009 Das 1.45 1.74 
12.5 2:15 1.10 127) 
SO 1230 0.80 1775 
5 aS 0.65 1.80 
20.0 0.85 0.40 1.90 
2255 0.60 0.30 2807 
25.0 0.30 Une 27 
205 0.25 0.22 2220 
30.0 022 QNOD 20 


3. Changes and Justification 


In order to achieve the desired pitch 


Home moments as listed іп table 


was not intended to modify the modeling of the 


components like seals, 


or fans, which surely would 


comparison with program 3, 


30 


sidewalls, 


angles one or 
1 һай to be changed. As it 
individual 
rudder, propulsion system 


endanger a qualitative 


the only moments left were rNWAVZ 
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STEED STATE BITCH ANGLES 
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PITCH ANGLE (DEGREES) 


x = IDEAL“ CURVE 
ЕБЕ ieee FURIES CURVE 
O = SIMULATION CURVE 
(PROGRAM 2) 
2.5 + 
x 0 
0 
>: 0 
p. F 0 
0 
0 
0 0 
0 z 
S t 
4 
= 
+ 
20 + 
t 7 
+ 
>< 
0.5 + 
+ 
+ Ж 
+ x 
+ ж 
%----%--------- %--------- %--------- %--------- %--- 
10.0 1550 20.0 250 30.0 


VELOCITY (KNITS) 
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and FMBUP assuming the aerodynamic forces found in wind 


tunnel test being correct. 


In figure 3 the correspondong forces and their points 
of attack are shown togetner with the according equations as 
used in the subroutines INCON and RHS. Following the 
calculations it can be seen that the monent FWAVZ is 
generated by the bubble wave making drag FXPWAV, which 
basicly is a function of the craft's velocity and weight. 
PEBAR can be considered constant similar to the pressure 
values in table III. in order to reach smaller pitch angies 
at higher speed FXPWAV would have to be increased 
substantially with the result that the water slope 
Correction term wATSLP and the thrust of the craft attained 
unrealistic values at lower velocity where the pitch angle 
tends to become negative as well as unforseeable changes in 
the dynamic response may be introduced. Therefore this 


possibility was rejected. 


ÀS nentioned before the plenum pressure can be 
considered constant. 30 the moment left Гог а possible 
pitch correction, FMBUB, could be used in this sense only Бу 
changing its lever, XCP, in a speed dependent manner. The 
method to determine the neccessary shift of the center of 
pressure (C.P.), in general, followed that outlined in Ref. 
1, the resulting values are given in table V and graphed in 
figure 4. During this procedure it was tried to match the 
pitch angles of the "three turns" curve. It can readily be 
Seen that there is an aluost linear relationship between the 
location cf C.P. and the craft's velocity. The following 
linear function used to approximate these data points was 
then included in the subroutines RHS and WAVES : 

XCPU = XCP - 0.034*Ux0.5921 + 0.3 
This new variable XCPU was used in the force and moment 


calculations. 
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FIGURE 3 
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FIGURE 4 
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Boo ee Onno pC enter Or Pressure 

(PEGQEAR 2, L£LESt change) 
Speed 00 
(knots) (reet) 
ПОО 10.36 
|255 10,25 
15.0 10.14 
(5 10.06 
2020 9.95 
2295 SEN 
25:0 92091 
یئ‎ ۹5 9.75 
30.0 5.68 


The updated steady state values are contained in table 
VI. Whereas the pitch angles approached the destinated 
values in a satisfactory way, the greater draft (апа 
consequently increased thrust) especially at higher speec 
may be regarded a setback. But as at this point it could ре 
seen that a similar correction would be neccessary for 
program 3 its impact on this comparison was not considered 


Crucial. 


Second Change: 

When caiculatlng the moment due to the bubble wave 
making drag, FMBUB, program 2 contains a correction term for 
I mno subroutine RHS : 

FHBUB = ABPB*(XCP-THETA*Z) 

As Z is used as a negative quantity (computed in subroutine 
INCON : 2Z= -ZS*DSO/12.0 ), FMBUB would increase uith larger 
pitch angles provided ABPB does not vary much. Under the 


small angie assumption for the pitch anale THETA this lever 








change generally is true, but in the above situation does 


not represent the position of C.P. of the XR-3. Presumably 
the location of C.P. on the 100B craft was almost underneath 
ЕНЕ септег оғ агауітсу (C.G.), with ZCP much larger than XCP 
and nearly equal to 2. In the case of the XR-3 for the given 
data ZCP is as large as XCP and considerably 


only twice 


smaller than the average 2. 


Table VI 


— —  -— — — سے‎ me * < s “a < m. s s + = шы - 


(Program 2, after first change) 


Speed Pitch Draft Pienun Thrust Rudder 
angle pressure 

(knots) (deg) 70 71٦ (1b) (deg) 
10.0 Ius DISSI 24.34 ПОБ T00 
12.5 do 1 7 8.18 20.46 855756 020 
15.0 0.94 8.10 PN TES 3500520000 
17.5 0.76 7.99 24.46 OS OO 
20.9 0.61 7.86 24.46 WB DE 
0025 0.48 ПЕЛА 24.46 486.60 0.0 
25.0 0237 70 24.46 SES OO 
2.5 0.27 2 7 24.46 0 9۹۹۷ ۵ 
30.0 0.19 7.18 24.46 715.64 0.0 

Therefore this tern was removed and by the foliowing 


function replaced: 

ХОРОВЕ ИЕП ХЕРИ, ZCP, THETA, PI) 
This corrected x-coordinate was than reintroduced in the 
subroutines WAVES and RHS, the letter using the new value to 
calculate the moment: 

FMBUB = ABPB*¥XCPC 
NUErunctvon SHXYAX = Shift of x-coordinate due to rotation 
around the y-axis - is kept quite general in order to enable 


its use for other cases. In this subprogram the radius of 
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the path that discribes the movement of C.P. is calculated 
and its angle with the X-axis. 
| HUNETION Z, ANGYAX, PY) 
Н = SQRT (Х**2+7**2) 
7'۶۹۹ 9+ GO TO 1 
NRE = Z7% 
ANGOLD = ATAN (ARG) 
The IF-statement avoids the calculation of ARG with a zero 
denominator. Depending on the sign of ANGOLD the true angle 
between H and X-axis 1S computed and the (pitch) angle 
ANGYAX suktracted: 
ET (ANGOLD СЕ 0-0) GOTO 2 
ANGNEW = ANGOLD + РІ -ANGYAX 
GO TO 3 
1 ANGNEW = PI/2.0 - ANGYAX 
GO TO 3 
2 ANGNEW = ANGOLD - ANGYAX 
mmemeerrected x-coordinate is finaliy found ру: 
3 SHXYAX = H*COS (ANGNEW) 


After introduction of this change the simulation 
program was re-run and the final results recorded іп teole 
VII and figure 5. Though this second change had oniy little 
influence on the steady state conditions as measured it was 
felt that it should stay in the simulation program because 
of its possible effect on the dynamic response where larger 


pitch angles may occur. 


The new moments are given in table VIII. They were 
attained under the same run conditions near steady state as 
those of table I. Due to the now negative moment FMBYB (а 
total shift by about 2000 lb-ft) and smaller pitch angle the 
bow seal increased its moment considerably. whereas the 
sidewalls reduced their effect as generally could be 


expected. 
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FIGURE 5 
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VELCCITY (KNOTS) 





(РоОодгаш 7, 


Table 


ーー ーー ーー ーーー Ow 


VEL 


ii ee ee 


after second change) 


Draft  Pienum Thrust Rudder 
pressure 

(in) (psf) (1b) (deg) 
8.31 24.36 405.74 0.0 
8.17 24.45 358.42 

8. 10 24.46 HUND O 
7.99 24.46 383.56 

7.85 24.46 428.20 

7.71 24.46 486.38 

ЕЛ 24.46 555.02 

7.37 24.46 631.98 

7.18 24.46 AG و‎ 


Table VIII 


simiy ol iom S (pitch) 


Speed Pitch 
angie 
(Knots) (deg) 
0.0 1.50 
5 1.21 
1220 0.98 
ШІ 5 0.79 
20.0 0.63 
22,5 0.49 
25.0 ШЕСІ 
27.5 (000217 
30.0 0:20 
(Program 
Conponent 
Name 
DS 
PSS 
FHS 
FMRUD 
FMP 
FMAED 
FMBUB 
FWAVZ 
` FMWAV 
Summary 


2, alter second Change) 
POSE Negative 
Joment fonent 
(EDS) (Bon) 
1227 E 
zc 61 
ج‎ po 
=== 83 
1792 UE 

650 — 
ーー ニニ 109 
— 277 

0 — 
3169 3169 
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the with the 


Puch 


further 


angle, 


For comparison program 3 


variables z-displacement, pienum pressure and 


thrust (starkoard) versus time were plotted 
9) 
JT. 


conditions 


(Plrotszszthroceh 
the initial conditions as given in table 
the the 


arting from 


enough in this part of state 
of 
conclusions concerning the dynamic behaviour already can 


plot 5 


study steady 


were primary interest, some general 


be 


drawn. additionally indicates a more damped pitch 


response due to the first change. 


Table IX 


— —— ーー= ーーーーー ーー ہے ہے‎ 


no en um “= 


(Program 3, original version) 


Speed Pitch Draft Plenun Thrust Rudder 
angle pressure 
(knots) (deg) (in) (psf) (1b) (deg) 
10.0 0.94 6.68 24.88 521.9 : 
12,5 0.68 6.48 DI 3 428.0 š 
15.0 0.47 622 aus? 392.20 
005 0228 SEO 24.81 38950 
20.0 0.11 565 2 ٦ 406.3 
22.5 0-69 O 24.86 449.3 s 
DOO 0898 j 6.03 24.86 492.7 0 
27.5 03:56 ٦ 24.87 543.6 
230.0 0.48 5.46 24.98 و‎ 71 5 
Hot o 22 - 10,06 feet 
В. РКОСКАМ 3 
1. Initial Runs and Evaluation 
The basic craft data used in Ref. 3 fortunately did 
not differ from those of Ref. 2 exept that the center of 
gravity (C.G) was positioned slightly more forward. The 
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results of the first simulation runs are summarized in 
tables IX and X under application of the same initial 
conditions, the pitch angles are graphed in figure 6. It was 
immediately apparent that in this program version, too, a 
correcticn of the pitch angle was neccessary before a first 


Comparison could be made. 


Steady State Pitch Angles 


に ーー ーー ーー 


(Program 3, original version, and Ref. 5) 


speed "ideal" X "Three Turns" Program 3 
(knots) (deg) (deg) (deq) 
10.0 2.45 1.45 0.94 
125 2.15 1.10 0.68 
0550 ШЕ) 0.80 0.47 
17.5 1225 0.65 0.28 
20.0 0.85 0.40 0.11 
22.5 0.60 0.30 0.69 
25.0 0.30 0.25 0.63 
21.5 0.25 ODD 0.56 
30.0 7522. 0.20 0.48 


2. Changes and Justification 


First Change: 


А deeper look into program 3 revealed that in 
Subroutine RES the speed dependent shift of C.P. of Ref. 1 
was still in use, whereas subroutines INCON, BOWSL and RHS 
already contained water slope correction terms iike Ref. 2. 
EL Uc existing calculation of the shift of C.P. was not 
BN sctory (the cause for the discontinuity at 22 knots 
remained unknown), these terms were removed from subroutine 


MMS anda temporarily replaced by XCP. Additionally all water 
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slope correction terms were revised or introduced in 
Subroutines INCON, SIDEWL, RHS, BOWSL and STNSL as outlined 
in Ref. 2 ain ali details. Besides these changes the sign 
errors in subroutine PROP of Ref. 1 discovered by Ref. 2 
were removed. A few trial runs, however, indicated only a 


negligible effect on the steady state values. 


The results of the following Simulations are reflected 
in table XI and figure 7. Despite the substantial changes 
made in nef. 3 the overali tendency of the craft's pitch vas 
almost identicai with that of the original version of Ref. 2 
malé a further proor that the introduction of WATSLP 


alone was not sufficient to reach proper pitch angles. 


"s awer ہے سو‎ ee ee = ъ= 


ーー = > 


(Program 3, after first change) 


Speed Pitch Draft Plenum Thrust Rudder 
angle pressure 
(knots) (deg) (in) (psf) (1b) (deg) 
10.0 1.97 7.05 24.90 405.67 
2.5 1.82 7533 24.92 puo 
15.0 S 7.09 24.91 334.35 
05 1.75 7202 24.84 347.26 
20.0 1.78 6.93 201 Р 
22.5 1.85 6.76 24.61 415.74 
25.0 1.97 6.49 u 459.64 
27.5 2.12 6.13 24.45 504.41 
30.0 2.15 5.93 24.35 555.91 0.0 


Note: XS = 10.08 feet 


Furthermore the center of gravity was re-positioned in 
order to match the conditions used in Ref. 2 (see II.C). 


This seemed permissable as both Refs. 2 and 3 used the sane 
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FIGURE 7 
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са configuraticn and attained a shift of C.G. by adding 
some load to the craft. The new location mainly effected the 
pitch angle as indicated in table XII. The simulation curve 
of figure 7 was lifted up by an almost constant amount, but 
the general shape was not influenced. Therefore, the task 


Ned to introduce an additonal correction term. 


Table XII 


— — — سے سے‎ «тре ہے کے‎ . z s. سے سے میت‎ s у чш сс «ж-ш. u често чете. 


Шіеоагат 1; агсег first change and Shift of C.G.) 


Speed Pitch Draft Plenum Thrust Rudder 
angle pressure 

(knots) (deg) (in) (pst) (1b) (deg) 
10.0 2715 7. 03 24.87 404.73 0.0 
15.0 1.90 OT 24.87 332.99 

2050 1.92 6.97 24.64 SUME 
25.0 P 6.46 24.43 455.68 0.0 
30.0 20 6.10 24.31 558.27 0.0 


At this point the corresponding moments near steady 
State were recorded in table XIII which enabies a direct 
comparison with program 2 under similar conditions (see 
table 1). The reduced draft and plenum pressure as vell as 
the changes made in Ref. 3 can be considered to account for 
the different moments FMBS, FÁASS, FMSW and FiBUB. Again the 
same variables as in program 2 were output (plots 9 through 
120 


05 





Mami TITI 


ーー 


(Program 3, original version) 


Component Positive Negative 
Name Moment Moment 
(1b-in) (lb-in) 

FMBS 850 کے‎ 
FMSS --- 324 
FNSH --- 3694 
FMRUD ーーー 83 
PIR 1049 --- 
FMAED 651 en 
ЕМВОВ 1831 mn 
FWAVZ =т= 280 
FMWAV 0 ہے‎ 
Summary 4381 4381 


Becond Change: 


The following modification was carried out the same way 
BEEN Nn -U*i1rst one of programa 2 (see III.A.3). The resulting 
values are represented in table XIV and graphed in figure 8. 
D педаг curve fit did not suffice a satisfactory 
approximation any longer a quadratic function was developed 
IND па1 form: | 

XCPU=XCP+0.001975% (U*0.5921-30.0) **2-0.97u4 

This new variable was used the same way to calculate force 
and moment. The resulting pitch angles besides the other 


Steady state values are shown in table XV. 
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a ہے ر‎ w Qj سے ہے سے چکھے‎ ш.ш» 


(Program 3, second change) 


Speed ХСРО 
(knots) (feet) 
10.0 10. 22 
1275 10.05 
15.0 85 
По ئ٦‎ 
20.0 9.56 
22.5 9-50 
250 9.46 
205 ӘЛЕ 2 
30.0 9.42 
Table XV 


ーー ーー Gm. ہے ہو‎ ーー ーー ere чы з ee we cee  — ーーーー ーー ーー ーー 


(Program 3, after second change) 


Speed Pitch Draft Plenun Thrust Rudder 
angle pressure 

(knots) (deg) (in) (psf) (1b) (deg) 
10.0 1.36 7.11 20.89 407.85 
12.5 1.03 6.89 24.86 348.81 
15.0 0.78 6.61 24.86 334.54 
7-5 0.59 6.31 24.85 345.26 

20.0 0.44 6.05 24.84 DAS 0 
225 0.33 5.81 24.84 410.41 

25.0 0.26 5.61 24.84 457.48 

27.5 0.23 5.44 24.85 511 GI 9٦ 
30.0 0.20 DIY 24.97 568.62 0.0 
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Third Change: 
E las alteration reflected the shift of C.P. due to 
E CIait's motion “as done by introducing the 


corresponding function 


ИТ. А.З. 
ШЕСІ pitch 


angles 


SHXYAX which 
table XVI and figure 9 show the 


this was only marginal a new curve fit 


once 


knowing the 


306 tene 


relationship it 


is described in detail 


final results 


slightly off the anticipated values. As 


Was not done, Dut 


would реа 


relatively short and straight forward procedure (see V). 


lable XVI 
Steady State Conditions 
(Program 3, after third change) 
Speed Pitch Draft Plenum Thrust Rudder 
angie pressure 
(knots) (deg) (in) (psf) (1b) (deg) 
10.0 1.38 7.16 20.81 407.95 
1185 05 6.93 24.84 384.90 
15.0 0.80 6.64 24.84 334.64 
in. 5 0.60 6.33 24.84 +4 ۹ ۹ ص‎ 
20.0 0.45 6.06 24.84 372.36 
2042 0.34 5.82 20.88 410.52 
P50 0.28 5.62 24.84 457.54 
275 0.24 с 24.84 512.04 0: 
ESO 0.24 SE 24.84 568.60 0.0 
Some the Steady state moments have changed 


XVII Due. Sto The 
moving C.P. most of the compensating moment is 


thats, TER 


considerably as can be seen from  tabie 
further aft 
picked up by the bow seal. It must be notified 
total moments of the corrected program 3 were 


than 


general, the 


с шошогасап Пу larger those of program 2 after 
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modification and would have had greater ci ect on the shear 


forces and bending moments which are not considered in 


this 


simulation(see V). The listed variables were also plotted 


for comparison studies (plots 13 through 16). 


Table XVII 


Summary of Moments (pitch) 


(Progran 3, after third change) 


Component Positive Negative 
Nane Moment Moment 
(lb-in) (lb-in) 
FMBS 2985 
IS SS == 
ЗЕК --- 2060 
FMRUD === 
ЕМР 1138 
FMAED 651 
FABUB BI 2094 
FWAVZ === 
FMWAV 0 
Summary 4774 4774 


۰ئ۱٦‎ ۸٠٠٥٢٠۱٢ AND RESULTS 


Although the data so far obtained allowed 


and conclusions additıonal information was 


neccessary 


sone comparison 


to 


demonstrate the basic differences between the two program 


versions. In table XVIII the interesting 


variables 


aere 


listed which were obtained in typical runs at 20 knots after 


all changes mentioned had been introduced. 


tine 


histories for pitch angle, draft, plenum pressure and thrust 


(starboard) give a good basis for a qualitative 


evaluation. 


Starting iron these trajectories a significant difference 
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can be seen in the pitch angle oscillations and the damping 
factor which will be investigated more thoroughly in the 
following chapters. The higher thrust in program 2 is 
certainly due to the greater draft so that under these 


consideration both programs reached reasonable steady state 


Conditions. 
7 ٤ 
Steady State Performance Data 
(Modified programs 2 and 3) 
Variabpies Units Program 2 Program 3 
Speed kt 20.0 ОО 
Pitch angle deg 0.63 0.45 
Time constant sec ЗО 3.06 
Frequency of 
oscillation Hz 0.48 0.82 
Draft in Lee 6.06 
Thrust lb 429220 377.6 
Air flow rates | 
Bow fans cu ft/sec 120283 20253 
Main {ans cu ft/sec 12122252 220 
Stern fans cu ft/sec ES 1279 
total cu ft/sec SS 8 5 2 
Pan power required hp 17.61 PRESS 
Actual fan power hp == 893 
Plenum pressure psf 24.46 24.84 
Leakage flow rates 
Bow seal cu ft/sec 172 10.4 
Sidewalls cu ft/sec 0.0 0.0 
Stern seal cu ft/sec 397 90 44.8 
voran cu ft/sec JOSE Бра 
Plenum volume CE EE е 399545 


Nevertheless, program 3 represents more realistic air 


flow and leakage rates in connection with the available 
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plenum volume as indicated in table XVIII (The difference in 
the voluaina basicly i: a consequence of different draft and 
pitch angles, too). The changes introduced had neglicible 
effect on the computation time in both cases, but program 3 
needed slightly less CPU time due to its better damping 


factor. 
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ПИТУ ОШЕА НОЕ 


(Dynamic Response) 
A. INITIAL RUNS AND EVALUATION 


In order to study the dynamic behaviour of the simulation 
models several possibilities were given to introduce the 
neccessary disturbance. Regular or irregular sea states in 
tabular form or by superimposed sine waves could have been 
entered using subroutine WAVES. Another way was to change 
speed by means of a thrust map of data points converted into 
a time history in Subroutine PROP. Finally rudder angles 
could have been applied following а straight run or by 
varying angles given as data points, too. The choice was 
determined primarily Dy the idea that besides the 
qualitative comparison of mathematical models a future 
validation of the results obtained should be made possible. 
Therefore, an easy~to-measure calm water situation was 
prefered. AS additionally new changes effecting the roll 
damping were to be incorporated the excitation through the 
speed did not seem very promising. Once the steady state 
conditions being round and used for initialisation a step 
rudder input (or a steep ramp input to come closer to tne 
real case) seemed the best solution offering a fine control 


of the disturbance magnitude. 


At a speed of 20 knots -actually the corresponding l 
initial thrust was kept constant - a rudder angle of 35 
degrees was found to give easily observable pitch and roil 
angle excursions for both models whereas larger values 
eventually caused water contact with the top cf the plenum 
chamber and stopped the program execution (controlled by 
Subroutine SIDEWL). 


Tne first Simulated turns under these conditions were 


recorded in plots 17 through 24 FOE both prograas 
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respectively (pitch and roll rates were changed into output 
variables in subroutine COLFIL and graphed in order to ease 
ПІ есъзша в 1 сп of damping factor ratios). As previously 
learned program 2 Showed a less damped pitch angle, but its 
roll angle response - after removal of the Superimposed low 
frequency component of undetermined origin = is slightly 
more damped than that of program 3. The speed dependent 
shift of XCPO therein may also be the cause for an increased 


tendency towards self contained pitch angle oscillation. 
٠۰۰٠٦۰٠ ۸۷ ۲2۸111080 OF SUSROUTINE SIDSEKNL 


Before the next changes are explained a Short 
introduction into the subroutine program SIDEWL is presented 
in order to give a better understanding of the possibilities 


and limitaticns of these alterations. 


The hydrodynamic and hydrostatic forces and moments 
acting on the sidewalls are found by appiication of slender 
body theory. This theory 1S based on motion of a foil in a 
viscous medium. The methods by which the lift and drag 
forces associated with moving underwater bodies are founda 
mM INMarily linear in nature (Ref. 4 and 8}, but heave, 
pitch, and roll motions cause the most Significiant 
nonlinear variations. Additionally there are certain 
nonlinear terms due to cross-flow drag, vhich are important 
for the case of very low aspect ratio lifting surfaces as in 
the case of the XR-3 sidewalls. The sidewall forces will 
then be composed of terms due to sidewall buoyancy anda 
Slender body hydrodynamic reactions and the effect of 
ii lOW drag: 


= PF ЕСЕ 
sidewalls slender body cross-flow drag 


In subroutine SIDEWL first the draft of or the gap 
beneath the sidewalls is calculated. The draft/gap is found 
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for each X-direction section after it has been corrected for 
craft roll and pitch angles and waves including water slope 
effect. In the presence of gaps the total sidewall leakage 
area ALSW is computed by summing the products of gap heights 


and sidewall section lengths for both sidewalls. 


In the following force and moment calculation the 
cross-flow drag terms and slender body theory terms are 
computed in Separate sections and finally summed to give the 
total values. In the case of the roll moment FK this 
calculation is of the forn: 

ЕК = (FZH(2)-FZH(1))*YSW-ZS*FY+FKD 
The slender body theory teras FZH for each sidewall (1=-port 
side, 2=starboard) are given by: 
| 25 —6Gm5C0-(U2UXxA03S#THETA-U*A33S*W+0*U 
Bes, pees UFAIISHPRYLSU 


where the expressions used are defined as follows: 


БСО = mass of displaced water by the jth sidewall elerent 
A33S = vertical added mass of jth sidewall element 

AD -XS 

BC2 = total sidewall vertical added mass 


YLSW 


N 


lateral distance of jth sidewall element from the 
centerline 
These components are either calculated in the subroutine 
itself or looked up in the sidewall tables computed in 
subroutine SIDETAB. Further details concerning the set-up 
of these tables are described in Ref. 4, the resulting 
values are stored in COMMON-block WAVTAB. The total lateral 
force FY is defined by: 
FY = FYH(1)+FYE(2)-FYD 
ТЕ represents the hydrodynamic lateral forces plus the 
corresponding cross-flow drag component derived as follows: 
ИПИЕ Ад OKU ま CVT XSSIR= ZS*P) 
where A22S is the lateral added mass at the stern given by: 
A22S = (RHO*.4*PI*DSS**2) /2 


with DSS beeing the computed sidewall draft. The final term 
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FKD is entirely due to cross-flow drag on both sidewalls 
each of them defined by: | 
ЕКО = FYD* (ZS-DSWAV/2) 

where FYD is the lateral force due to cross-flow drag and 
DSWAV is the draft of the jth sidewall section corrected for 
waves. In the computation of FK the lateral force FY is 
used without the drag component which is addea later. In 
general it must be said that all forces more or less 
influence each other so that a change in the roll 
characteristics neccessarily must alter the pitch and уан 
behaviour or vice versa as could be seen when intrducing the 
pitch angle correcticns or in the first turn runs where the 


pitch angles increased significantly. 
C. CHANGES AND JUSTIFICATION 


The following changes were recommended by Bentson [Ref. 
ENDE ev consisted of two major blocks and a minor addition 
carried out in subroutine SIDEWL. The letter dealt with the 
vertical added mass at tne stern, A22S, which is recomputed 
for the case of negative pitch angles using the corrected 
draft at the bow, DRBON: 

DRBOWZDSS- (XX (J,Nt1)-XSS) * THETA 

Ir (DRBOW.LT.0.0) DRBOW-0.0 

A22S= (RHO*.8*PI*DSS**2) 72. 

Poche thet ta.) AZZ5—.6*hHO*PI*DRBOX=DRB OW, 2. 
XX is the longitudinal distance of the nth sidewaii station 
ШЕСІ (пе boy. The factor .8 in the calculation of A22S was 
changed to .4 in order to remain consistent with the ھ٥1‎ 
given based on the definitions derived in Ref. 4. The 
effect of this change, however, was not observable because 
the Simulated pitch angles were most of the time positive 


except a few small negative transient peaks. 


The next major change was a correction of the vertical 


force FZH. As the outside edge of the submerged porticn of 
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the sidewall is not vertical near the bow an additional 
upright force component lis generated due to the  deadrise 
٣۰۰ ۰ ٢۰٢٥ب‎ of the lateral force FYH, in algebraic terns 
defined by the product of the lateral force and the 
cotangent of the deadrise angle (average angle between tne 
horizontal and the inclined outside edge of the mentioned 
sidewall portion). In FORTRAN IV Language the computation 
proceeds as follows: 

CTNDR=0.0 

ШЕ (робите (0) CO TO 22 

CTNDR= (BS-BB(1))/DSS 

Api; 70) C NbR=0.39391‏ ۱۱ت 

22 CONTINUE 

FZH (39) =FZE (9) +PM1*FYH (9) ¥CTNDX 
BS and BB(1) are the lower and uvper beam widths of the 
submerged sidewali portion, CTNDR stands for cotangent, and 
PM1 is a unity sign factor depending on the sidewall under 


duseWesslon. 


The last addition had tne purpose to improve the roll 
damping of the model. It vas developed from the idea that 
when the craft is rolling the sidewails generate asymmetric 
waves depending also on speed which take energy out of the 
ship and, therefore, have a damping effect on the roll 
motion. Without discussing each single step of this rather 
involved change the general procedure was to evaluate hew 
E bers corrected for these waves and from there to 
find correction terms for the vertical added sidewall mass 
BC2. The previously determined roll moment FK was finally 
adjusted by subtracting the resulting counter moments 
including an experimental factor. The complete change was 
done as follows: 

Mio — Ate و‎ 3۶ 

ZONI (SLGN( Is, DSS) +12) 7/2. 
170 1 

DS=Z+ZS 
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STDS CARE XS) THETA 
ID=1.+(DSR*12.-SDS) /DDS 
ID=MAXO (MINO(ID,NDS),1) 
DDSR= (1D- 1) *DDS+SDS 
ID1-MINO (ID*1,NDS) 
DID=(DSR*12.-DDSR) /DDS 
BOO NCO سس۶‎ ٦ 
BC2-BC24DID*(AC2(1,ID1,ID)-BC2) *DIP* (AC2(1,ID,IP1) 
-BC2*DID* (AC2(1,ID1,IP1)-àAC2(1,1D,IP1) 
-AC2(1,1D1,1P)+BC2)) 
BECERT-UOSGUTSWRTYSHdTBCO*P/PI 
For the shear and bending moment calculetions not contained 
in this program version the vertical force components of the 
individual sidewalls were also corrected: 
FZH (1) =FZH (1) *8.*YSW*B8C2X*P/PI 
FZH (2) =FZH (2) -8.*YSW*BC2*P/P1 


Reference 9 finally indicated the option to inciude the 
calculation of the spray drag of the sidewalls when 
Operating at high speed with minimum draft. Whereas this 
function originally was developed for the  100B craft with 
considerably smaller draft values down to one inch at about 
80 knots the XR-3 operates with five to six inches or nore 
and, therefore, the effect of the spray drag vas considered 


negligible. 


VS COMPARISON AND RESULTS 


Starting with progran 2 the new roli and pitch 
characteristica were recorded in plots 25 through 28. 
Besides a reduced overshoot of the pitch angle by about 20 
per cent the roll angles reached only about half of their 
original peak values while the non-linearities of tne 
sidewall changes became quite apparent. Therefore, the time 
constants could not but be estimated and seemed to have 


increased slightly for both motions. Although the overall 


9 











response could be considered an improvement "an order of 
magnitude difference in the roll damping value" аз pointed 
out in Ref. 9 and meant for the Bell 100-ton test craft 


could not be experienced for this model. 


For program 3, whose time histories were graphed in plots 
29 through 32 the changes introduced were not as well 
favourable. The roll angle overshoot reduced only by 209 per 
cent whereas an equally diminshed pitch angle Seemed 
considerably less damped and showed the increasing trend for 
self-contained oscillation as indicated earlier. However, it 
must be noticed here that program 3 originally vas not quite 
ав underdamped as program 2 and, therefore, could not be 


effected by these changes by the same amount. 


In order to attain an optimal model based on tne previous 
results and considerations the two major changes made in 
subroutine  SIDEWL and their individual inpact on the 
simulated performance were studied more thoroughly afte 
ili a few modifications. In the calculation o£ the 
uucead force due to deadrise projection of FYH the 
Statement 

FZH(J)-FZH(J)*PMT1X*FYH (J) *CTNDR 
was replaced by the following expressions: 

FZHOLD (J) -FZH (J) 

FZHDRP (J) =PNI*FYA (J) *CTNDR*PRO4OI 

FZH (J) =FZH (J) +FZHDRP (J) 
These statements allowed to print the individual force 
components if desired and to vary the component value due to 
deadrise projection by means of the program card deck using 
one of the modification settings (see Appendix А). Similary 
the roll damping computation was re-organized: 

FROLD=FK 

ET pi en Oi ZY oWKY sre BCZ*EPs7 PI 

FZH(1)=FZH (1) *PROMO2/2. *YSW*BC2*P/PI 

FZH (2) =FZH (2) -PROMO2/2. *YSW*BCZ*P/PI 


60 








Finally a  WnITE-statement was added using PROMO3 as print 
Switch for VAL(1), ЕЛ ДОБР (Ту» FZHOLD (2), FZHDEP (1), 
FZHDRP (2), FKOLD, and FK. The influence of various settings 
of PROMOI and PROMO2 - here simply called deadrise factor 
and roll damping factor - could be studied on the following 
plots 33 through 60, the meaning of the R.D. number is 


ЕИшегі?есй in table XIX. 


Tanie zix 


Listing of Damping Coefficients 


جح — > — 


(Program 3) 


EDS number Deadrise Factor ۲۰٠۰٠۰٠۷۹ ٠٦۷٠۷٠٠٠٦۶ تق‎ ٥٣ 
1 ie 10% 
2 0 16. 
3 т: 22. 
4 3; То. 
5 (Є 48. 
6 Ü 1. 
7 Io си 
8 5 а 


Using program 3 in these simulation runs it could be 
shown that the response was improved by lowering the 
deadrise factor and increasing the roll danping factor as 
Di] plots 57 through 60 under R.D. number 8. These 
values were chosen to indicate an approach method only, more 
extended Simulations with a finer factor spacing could 
possibly result in an even better damped behaviour. This 
step, however, should be done in conjunction with a 
validation of this model in comparison with measured test 


data. 
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VEC ONCU CLONS AND. RECONBENDATIONS 


a = < s 


A great deal of time - more than anticipated at tne 
onset - was required for familiarization with the computer 
models and for the implementation of comprehensive and 
mnossonable initial trim conditions. For future studies 
Appendix A and the example program printout as well as the 
presentation of many tabulated data may help to reduce this 


introductory phase effectively. 


The simulated roll and pitch motions, in general, vere 
improved aithough some of the changes made need further 
investigation in connection with validation tests which 
should be performed with program 3 as the presently nore 
advanced version. In case of future changes effecting 
longitudinal and athwartships stability it is recommended 
that the shift of center of pressure is verified as outlined 


in Appendix A before concluding tests are undertaken. 


These further Studies may concentrate on the 
investigation of the moment due to the product inertia term 
AIXZ which has some importance in maintaining stability in 
roll and pitch for larger rudder angles. The use of the 
discrete mass distribution option (Block 2) can be helpful 
for this case. 

Anotner field of interest can be the air pressure 
distribution in plenum chamber and both seals and its 


ШІ Тепсе оп transient motions. 


Perhaps the most significant improvements may arise from 
a comprehensive study of the modelling of the sidewalls. 
Besides the calculation of the vertical force due to 
deadrise prcjection of the lateral force - the sidewall 


wetted surface shape of the Bell 100-ton craft is almost 
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uniform from bow to stern, whereas the XR-3 is more formed 
like а wedge widening at the stern - the impact of flat 
portions on the bottom of the sidewalls as indicated in Ref. 
3 should be investigated as they become effective planing 
surfaces at higher speed (20 - 25 knots) but are not 
considered in the simulation program so far. The need to 
lower the computed cotangent value indicates that the 
present calculation of the deadrise projection does not 
enough consider these structural differences. A pitch down 
due to a more negative sidewall moment would also reduce tne 
neccessary shift of the center of pressure and thus the 


destabilizing effect on the transient response. 
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INTRODUCTION 


This users manual is a short form of that given in 
reference (1) which should be referreä to in case of severe 
programming difficulty and for a deeper study of the 
susroutine prograns involved. The version given here 
reflects the changes iade in the meantime by Refs. 2 and 3 
and by the author of the thesis study the manual is attached 
to. Its objective is to introduce the reader to the general 
program flow, to the neccessary requirements to establish an 
own input deck for various types of craft runs and to the 


ШЕБІ options. 
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PROGRAM ORGANISATION 
Within a given run the program proceeds as follows (see 

Exunre 10): 

- First the main program calls the input and initialization 
subroutine INCON. The program then proceeds to calculate 
the motion time histories by calling the various modular 
subroutines апа the integration routine for as many time 
Bene as is neccessary to complete the run. During the 
Course of the calculation the various subroutines will 
output those variables selected by the print option 
switches at the appropriate input print interval. In 
addition, data may be written on Scratch files to be used 
for summarizing the output and doing bending moments and 
Shear calculations (at the present time the shears and 
moments subroutine has been removed). 

- When the run is completed, the program, depending on the 


options chosen, may then print and/or plot the output 





summary. After these steps are completed the main progranm 
sl to INCON to read data for the next case. 
- The above procedure is repeated until all cases are 


finished. 
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FIGURE 10 
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CONTINUED 
FROM Bro АСЕ 
x V 
RHS 4------- 4 BOWSL | 
(K1-K5) 4---- — — 4+ سب سب -“ ب‎ шь 
ا‎ STNSL | 


ғ mm am чии” «ие “Um ود‎ Um s s me 


= — ーー 


ж. «ки» «ки» ms a e e SEE Q M ow 


a am Homo o O та» MU “< e e (min e ана s 








- 3 ҒАМ [e EUM | 
| INTGRL | 
v 
| MAIN | 
v 
МООР. |4----- 3 PLOTP | 
ーー | DRAW | 
t 
E S m 
| 
END 
OR 
NEXT RUN 


* - [LIBRARY SUBPROGRAMS 
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ーー = に 


А5 гип оп ап IBM 360 computer under OS/360MVT version 
Ше саде 20.6), on W. R. Church Facility, NPS, the program 
requirements are as follows: 

TOTAL CORE 


w =m a mo e ۔‎ — + = u жу 


Program length Ре 


Core size requireed for execution nr 


The total core required for a given step includes the 
program, system routines, data storage, etc. The time 
required to run a single case depends on the nature ог the 
Case and the length of time to be simulated. The ratio of 
required execution time to real time has been found to be as 
low as 1/3 for some calm water cases. Higher ratios occur 
where the transient curves are rather sharp as during tne 
turn maneuvers in this study where 20 seconds run tine 
needed about six minutes using FORTRAN G compiler. However, 
Mitchell [ Ref. 10) Zound an average improvement of 


thirty-seven per cent CPU time by compilation in FORTRAN  H. 


The program currently uses the following FORTRAN Data 


Set Reference Numbers (Unit Numbers): 


Unit Number Usage 
1 Beraten tile Тос COLYIL plotting package. 
2 ratch file for Lateral plane output 
3 Summary 
5 ШЕШ ШИЕ 
6 Output 


10-15 Sidewall Integral Taole data 

IBMF Intermediate file containing time histories 
required by bending moment calculations. The 
unit number IBMF is read in using the 01501 
input card not contained in the present 


program version. 
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INPUT 

The input logic is subdivided into twenty-two "blocks". 
ШЕСИ block being for a particular subroutine or function 
EE Stern Seal, sidewall, etc. Switching to the proper 
E Lock іс done by means or a control tag on the input 
cards. Once control has been transferred to a particular 
block additional input may be read by specific format 
statements if necessary, however, the card which contains 
the control tag may also be used to input up to seven 
variables for the selected block. The format of the input 
card used for switching is described later, a sample data 
week. tor the last simulation run (R.D. number 8) is copied 


Eu end of this appendix. 


The опто tage Sed to transfer control to the 
appropriate block is of the form OZN where M and N are 
Ser Constants, i.e., if the control tag is 015 the 


program will go to statement 1500. 


The option tag is loaded into a variable called IOPT and 
allows different branches of logic within a given input 
block. If no branches are used within the input block the 


Option tag may be zero. 


The input variables are stored in an array called TEMP 
and may be set to a given parameter by using an arithnetic 
replacement statement in the appropriate block (see various 
mocks in attached listing for examples). Print switches 
bontrol the printing by each subroutine of its forces and 
moments. This option is very usefui in debugging the 
program. Normally it provides too much outout and is too 


cumbersome. 
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OUTPUT 
The output generally consists of the following parts: 

- The data input deck is first copied headed by the title 
ща ING OF INPUT DECK". Each block and option is printed 
as actually puncned onto the cards. These data are written 
in the main program prior to any test and, therefore, 
appear in ail runs. 

- After the data are read by subroutine INCON a more 
explanatory Summary is output entitled "SES MOTIONS AND 
LOADS PROGRAN® followed by the description written in 
Brock 18. This table contains most of the input data and 
first computation results like inertia matrix and xx- and 
E drrayvs oÍ sidewalls and seals. In cases of input errors 
a corresponding message follows and execution stops. 
Likewise all other problems preventing a proper completion 
LC current run are indicated at this place. 

EF plots are requested they are either printed іп the 
ہہ ہک‎ part ог, for CALCOMP plots, their titles and 
۰۰۰٦٦٣10 7 are Summarized using subroutine COLFIL. 

- In the final portion the tabulated variables as desired 
are summarized in one or two divisons followed by the 
message: 

CONPLETAD ALL RUNS 


On tne following page the possible variables for these 
plots or summaries are listed. A change to other variables 
needs only moderate program modifications іп subroutines RAS 


and COLFIL, eventually in some others, too. 
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— —  — > — — ` < а — سے‎ m= Q w rn m ہے‎ 


The names mentioned below are used in subroutine COLFIL, 


but may have different names in other Subroutine programs. 


Number 


01 
02 
03 
04 
05 
06 
07 
08 
09 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 


Nane 


TIME 
ETA 

2 

THETA 
PB 
ВОКАСС 
АСС 
FANPWR 
PHI 
BETAS 
ACCLAT 
U 
TRADIUS 
VOLP 

X 

Y 

QIN 
QOUT 
GFXXX 
ЕХРИАУ 
THSTS (1) 
THSTP (1) 
QDEG 
PDEG 
RDEG 
DELES 


موت فصت س m E‏ چ ہے کے چ کے تاس 


Independent variable time, seconds. 
Wave height, feet. 

Z displacement, inches. 

Pitch angle, degrees. 

Plenum pressure, psf. 

Bow acceleration, ft/sq sec. 
C.G. acceleration, ft/sq sec. 
Fan power, horsepower. 

Roll angle, degrees. 

Yaw angle, degrees. 

Lateral acceieration ft/sg sec. 
Speed through the water, knots. 
Turn radius, feet. 

plenum vojume, cu ft. 

X displacement, feet. 

Y displacement, feet. 

Air ilos rate in, cb o Lt/sec. 
Aur lo Na e опт me) ty GG 
Net force in x direction, lb. 
Wave force in x direction, 1b. 
Starboard thrüúst; Lb. 

port thpust, Ib. 

Pitch rate, deg/sec. 

Roll rate, deg/sec. 

Yaw rate, deg/sec. 


Rudder angle, degrees. 
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МІ РКОСНКАН АНПр SUBROUTINES 


Ñ — m "w a a an‏ سے سے 


Hain Program 


Discussion - The main program contains the logic for 
interconnecting the various  subroutines. In addition, it 
compares tne running value of time with the finish tine, 
calculates the next value of time for printing after each 
ШІП (іле, etc. The main program also contains the logic 
Ecc fixed coordinate trajectory (x, y, U ) by the 


Propazoidal rule. 


ШЕШЕ The output from the main program is controlled by 
the print option switch ITRAJ. If the value of tais switch 
is 1, the program will print the current values of the tine 
and the craft velocities and dispiacements (angular and 


translational). 


BLOCK DATA 


a a s s ee ee ы чь» Өл» ج کے‎ 


contain all COH4ON-statements used in various subroutine 
programs. All values therein are initialized with "0.0" 
sone the first data are read into the program. They keep 
their current values between consecutive runs (controlled by 
ПОО Л) as far as they are not recalculated due to a 
change of input values or within the subroutines themselves. 
Various constants, 200, are transferred by 
COUMON-statements. 


Subroutine INCON 


Discussion - Subroutine INCON contains the logic for the 


— — سم‎ omy ーーーーーーーー ーー ニー 


ЕЕВС Oof all input data, the initialization of variables 
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ana for the initiation of new cases. The subroutine is 
programmed so that the input data is non-structureä i.e. 
there is no set order of input cards from one read statement 
to another. However, a given read statement may have more 
than one input card. In runs with multiple cases, only those 
input parameters which are to be changed for the new case 


need be read in. 
Usage - CALL INCON 


utput - The output from INCON is discussed earlier. 


Subroutine DMINV 


Discussion ~ Subroutine DMINV is a matrix inversion package 
used for inverting the mass matrix. The method used is tne 
Standard  Gauss-Jordan technique. The determinant of the 
matrix 15 also calculated. If the determinant is zero the 


matrix is singular and the program will stop. 


Usage - CALL DMINV (A, N, D) 


= ーー oe ao 


A = The input matrix, destroyed in computation and replaced 


by resultant inverse 


N = The order of matrix (maximum = 6) p 
D = The resultant determinant of the matrix 
Subroutine NAVES 
Discussion - Subroutine WAVES calculates the wave forces and 


ーー p Zr mn | шш 


nls acting оп the craft. It also generates the wave 
heights at the various stations around the seals and 
sidewalls, as well as the bubble volume lost due to wave 


elevation. The sea state for irregular seas is computed in 
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M っ 


WAVES ky adding together a series of regular waves with an 


appropriate distribution of anplitude aná frequency. 


Usage - CALL WAVES 


ШЕШЕ - Тле output fron WAVES is controlled by print option 
switch IWAVES. When this switch is set equal to 1, WAVES 
Will Print the wave elevation around the seals and 
sidewalls, the wave elevation at the center of gravity and 
the volume reduction of the bubble plenum due to waves as 
well as the total forces and moments acting on the craft due 


to waves. 


Subroutine SIDEAL 


KT SIOR = Subroutine SIDEHL calculates the forces aná 
moments acting on the craft due to the sidewalls as well as 
the leakage flow rates associated with any gaps which open 


he sidevalls due either to craft motion or vaves. 


Me CALL SIDENL 


Ж РЕ - The output from SIDEWL is controlled by the print 
Ш ПОП Switch ISIDHL. When this switch is equal to 1, the 
program will print the sidewall gaps, immersion depths and 


total forces and moments. - 


Subroutine PROP 


Dascussion - Subroutine PROP is used to calculate forces and 
noments on the craft due to the propulsion systen. The 
Subroutine has the capability of uneven thrusting and 
receives this information by using the function program FG1 


to linearly interpolate a map input in block 16. 
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Wsage - CALL PROP 


ШОШ The output from PROP is controlled by print option 
ЖЕР UNPBROP. If this switch is 1, the program will print 


the calculated forces and moments. 


Subroutine RUDDER 


Discussion - SUSTEOULINGE “RUDDER contains the logic for 
calculating rudder forces and moments as well as rudder 
motions. Rudder motions are calculated by empioying the 
function program FGi to linearly interpolate an optional map 


input. 
Usaqe - CALL RUDDER 


Dutput - The output from rudder is controlled by print 


ION Switch IRUD. If this switch is 1, the program ‘twill 


print out tne total forces and moments due to the rudder . 


Subroutine AZROD 


Discussion -Subroutine AEROD calculates the aerodynamic 
forces and moments on the craft. The subroutine currently 
uses curve fits of the aerodynamic coefficients vs sideslip 


angle. 
Usage - Cali AEROD 
meme — The output from AEROD is controlled by print option 


Bohn  IREROD. ТЕ this switch is 1, the program will print 


the total forces and moments. 
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Subrourine INTGRL 


Discussion - Subroutine INTGRL is used to integrate a systen 
Of first order ordinary differential equations. It uses a 
variable time step technique based on the Runge-Kutta-Merson 
algorithm. А fixed time step Runge-Kutta algorithm is also 
available. The program will stop the calculation if tne 


time step becomes smaller than 10-9зес. 


Usage - CALL INTGRL (TIME) 


Description cf Parameters 


a a a e w ae‏ کد سد سے کد ل سد دد کے کد ٹک کک 


TIME = current value of the independent variable 


Subroutine RHS 


Discussion - Subroutine RES is the subprogram containing the 
FORTRAN expressions for the right hand side of the system of 
Meroe Order differential equations. It contains the logic to 
calculate the bubble volume, area and pressure as a function 
of time as well as the logic which sums the individual 
forces , moments and leakages from the various craft 
 Гпеп(гс (о get craft totals. RHS also contains the 
statements used for the writing of the scratch file for the 


۲٣۰٦٢٦۷٢۹۶7 output. = 


a 
Im 


Meme CALL RHS (VALUE) 
Discription of Parameters 

rue = The array containing the values of the right hand 
side of the differential equations as calcuiated in 


RHS, i.e., for a system of equations of this form 


ООШ (yA), t) the VALUE array IS used “to. store 


F 





the values of f(y(i),t). 


BEN put - The output from RHS is controlled by print option 
Ach IRHS. When this switch is 1, the subroutine will 
print all right hand side forces, moments, accelerations, 
pressure, air flow rates and fan power besides the complete 
VAL- and VALUE- arrays. This print switch is the most 
efficient one for debugging purposes but also computer time 


consuming. 


Subroutine BOWSL 


Be euss3on = Subroutine BOWSL calculates the forces and 
moments due to the bow seal and the leakage flow rate 
associated with any gaps which open up under various seal 


Pee tions. 


Usage - CALL BOWSL 


ИЕКТЕН Вон. ТЯ this switch is 1, the program will print 
the gaps and wetted lengths at the various stations on tne 
bowseal as weil as the total forces and moments dve to tne 


bowseal. 


Subroutine STNSL 


EM cussjon ~ Subroutine STNSL calculates the forces and 


m= w x s ہے‎ чш» с ーー ہے ہے‎ 


moments acting on the craft due to the stern seal as well as 
the leakage flow rates arising from any gaps which open 


under the seal. 


Ge = CALL STNSL 


© کے سے ーー‏ 
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AER. ISTNSL. When this switch is equal to 1, the program 
Will print the seal wetted length, and gaps as a function of 
Station as well as the craft forces and moments due to the 


Ern seal. 


Subroutine FAN 


Discussion - Subroutine FAN is used to calculate the inflow 
to the plenum, stern seal, and bowseal from the fans for a 
given pressure differential across the fans. the subroutine 
uses tabular data and allows for changes in fan speed as 


well as the number of fans. 


Usage - CALL FAN 


〇 


ШОО Т The output from FAN is controlled by print option 


| 


٠٠.٦ .5ظ‎ If this switch is 1, the program will print the 
flow rates and pressure differentials across the bowseal, 
cushion and stern seal fans, aS well as ideal fan power, 


actual fan power and fan efficiency. 


510700 0 68 COLFIL 


Discussion - Subroutine COLFIL is used to format the program 
٦٣۰٠۰. Tt Can provide up to 16 variables listed in two 
tabular summaries of 8 variables each. The Subroutine has a 
plotting package that permits use of either the print plot 
ca lcomp (continuous line) option. Multiple curves may be 


plotted and up to 10 graphs may be drawn on each run. 
Meade -CALL COLFIL 


Ment — The output from COLFIL is discussed earlier. 
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Subroutine SAM 


Discussion - Subroutine SAM is supposed to calculate tne 


shear and bending moments, but is not used in this program 


version and replaced by a dummy subroutine. 


puruction ٦۹٦ 


Discussion - The purpose of FGÍ1 is to perform the task of 
evaluating a tabuiar function y = f(x). the proceäure used 
is linear interpolation between data points. If the input 
value of tne independent variabie does not fall in the range 
of its input tabular values the appropriate table extreme is 
held. The function y must be single valued and the tabular 
values for the independent variable must be stored in 
Ascending order i.e. XTAB(J)XXTA3(U*1). The increnent 


Between tabular values need not be constant. 


Usage -Y = FG1 (X,N,XTAB, YTAB, IMEM) 


Description of Barameters 
Y = The output value o£ the tabular function 
X = The value of the independent variable to be used to 
evaluate the function y = £(x) 

N = The number of elements in the XTAB and YTAB arrays 
faye = The name of tne array that contains the values of 
the independent variable in ascending order 
IE Tie name of the array that contains the values of 

the dependent variable corresponding to the values 
of the independent variable stored in XTAB 
ТИЕМ = Dummy variable used to store location of last entry 


nto table, to be used as initial point of. search 


at next entry 
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Function SHXYAX 


EEScussrion - Function SHYYAX calculates the shift of the 
x-coordinate of a given location due to rotation about tne 
y-axis in the local reference frame. It is used to attain 


тесте levers the moment calculations. 


Usage - XC-SHXYAX(X,Z,ANGYAX,DI) 


| 


mere mes we ne سے ہے ہت‎ ee = ーー ーー cc سے لے‎ яж» <<< -— کے‎ 


€ = Corrected x~coordinate 
X zu rocg9)nal x-coordinate of given location 
2 amori ginal 2-coordinate of given location 
ANGYAX = Angle of rotation about y-axis 
рт = Constant pi 
۶۰۰٢۰۰۰٢ ons, Tiıand 22 
Discussion = Functions Ши and T2 are two small 


FUNCTION-subprograms used to calculate various trignomsetric 


relations used by subroutine WAVES. 
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The names mentioned in the following summary are used in 
subroutine INCON. Not all of them are transferred as read to 
other subprograms. An investigation of the corresponding 
COMMON-Statement may give further hints, but care must be 
taken as those names may be not identical in all 


subroutines. 


ーー ーー ーー ーー ーーーー ーー ーー ーー ーー ーー ーー ーー ーー ーー  — ہے سے سے‎ = 


Option tag 1 = Time steps 
Card Column Fornat Name Entry 
1 = 5 т 00101 
B 15 F10 STIER Start time (generelly zero), 
sec. 
16-25 F10 FTIME Finish time, sec. 
6-35 F10 DELO Initial integration time step, 
sec. 
6 5 F10 DELPNT Print interval, sec. 
46-55 F10 TPRINO print start time, sec. 
Optlon tag 2 = Print switches 
Card Column Format Name Entry 
1 1-5 I 00102 
2 5 I 08 ۳ Print switch for lateral 


acceleration and integration 
time step. 
10 т IVEL Not used at present. 


15 m ITRAJ Print switch for main program. 
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30 


39 


40 


45 


motjon tag 3 


Nane 


LSI DHL 


IDOWSL 


707 


IWAVES 


7 0 


I PROP 


IAEROD 


PRES 


Entry (cont) 


ры Ssywiechs for süubraurıne 
SIEDEML. 

Print switch for subroutine 
BOWSL. 

Print switch for Subroutine 
JINS: 

۶۲۰۰٠۰۰۰٠٣٠٠٠٠ ٠٠ ٠٠٠٠٦٦٦ ۶97 706 
WAVES. 

Би мен мек” "тог subroutine 
RUNDER. 

Print switch for subroutine 
EROP., 

Leite switche fon subroutine 
AEROD. 

Prime switch Тос subroutine 


RES and FAN. 


= Integrator tolerances 


Card Column Format 


— | | “se e — —— —  — — 


Í 0 
2 ШЕ 2 
3-4 


Name 


Entry 
00103 


Number of integrators to be 
used (maximum of 15). 
Controle tad ос» integration 


procedure. 


00 = Original Runge-Kutta 
procedure. 

01 = fixed step size 
Runge-Kutta procedure. 

02 = variable Step size 
Runge-Kutta procedure. Only 
the minimum step size is 
output. 
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ШГП column Format Name Entry (Cont) 


EU 1-80 8F10 TOL (J) Integrator error tolerance, 


one per integrator. 


ee часто 07س سیب‎ 


The usual procedure for selecting integrator tolerances is 
to choose values on the order of 10-5 to 10-8. For a new 
eraft the following procedure is used: 

Make a run with a given set of tolerances; halve the 
mererances and rerun. If the results do not change 
drastically, try doubling the original tolerances and rerun. 
As changes are noted in the integrator outputs (i.e. craft 
velocities, etc.), hold the tolerances for those integrators 
fixed and vary the others until the largest  tolerances 
possible are obtained. Once chosen for a given craft the 


tolerances rarely have to be changed. 


Option tag 4 - Tabular summary of output 
Card Column Format Name Entry 
1 15 I 00104 
2 5 I ТҮЗІ Print switch for summary one 
erint = nor print, 
10 1 ILATRL Print switch for summary two 
T eprint; 0 = no print: 
15 I NVD Not used at present. i 
20 I NVI Not used at present. 
25 ШЕ NLD Not used at present. 
30 1 NLI Not used at present. 


| Option tag 5 - Program option switches 








Card Column Format Name 


+ —— سس چیہ جت جه سے ج r‏ — — — 
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Card Colunn Format Mane Entry (cont) 


—— = — — سے سی ےہ ہے سے‎ ш єк 


1 ше 15 F 10 23 OR Switch TSF lateral plane 
motions. 
1 = only lateral plane motions 
0 = six degrees of freedcen. 
16-25 F10 TSRGE Switch for constant speed. 
1 


0 = surge velocity allowed to 


surge equation not used. 


vary. 
26-35 F10 TERIN EUCH ise t Pathe 
1 = thrust varied to give 
| proper speed. If this option 
| is selected, then there must 
| ре 3017 ІШ Block 16 7 
estimate of thrust and a 0.0 
for the number of data pcints. 
0 = thrust as input. 
E545 F10 PD A Switch for membrane study 


concerning plenum air bubble. 


Option tag 6 = Program modification settings 


Card Column Format Nane Entry 
1 Eo’ I 00106 
5-15 F10 PROMO1 Program modification setting 1 
(temporarily used as deadrise 
factor). 
06-5 F10 PROMO2 Program modification setting 2 
(temporarily used as TOLL 


damping factor). 

26735 F10 PROMO3 Program modification setting 3 
(temporarily used as print 
Switch for vertical forces and 
roll moments in subroutine 


SIDEWL) . 


85 





Card Column Format Name исту (cont) 


1 82.145 F10 PROHCU Program modification setting 4 
86-55 F10 РКОМО5 Program modification setting 5 
0-65 F10 PROLOG Program modification setting 6 
05 F10 PROHO Program modification setting 7 


Block 2 Hass and inertial Input 
Option tag 1 =1 Summary mass properties card 


Note: In this option shears and moments will not pe 


مسوم وس اس < — 


Card Column Format Name Entry 
1 1-5 I 00201 
6-15 F10 WEIGHT Total eral Welt, pounds. 
16-25 F10 5 Longitudinal center of 
gravity, feet forwara of 
transom. 
26-35 F10 ZS ۷۰۰۲۰۲۲۰۷٣۹ ٣٢۴٣٣ آ6 7ت‎ 0۶۰۲۴۷٣٦٠۷ 


feet above baseline (keel). 

36-45 F10 NY Mass noment of inertia about x 
AKL, oclug- sd ft. 

46-55 F10 түү Mass moment of inertia about у 
axic slug so ft. 

56-65 F10 А122 Mass moment of inertia about-z 
axis, SLlug~-sgqg rt. 

05-75 EIO А1Х2 Mass moment of inertia about 


AIS Sl] SG ft: 


| Option tag 2 = Discrete mass distribution cards 
Card Column Format Nane Entry 
1 1-5 I 00202 
ENS. 1-10 F10 AMI (I) Discrete weight, pounds. 
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Bard Column Format Name CATE CONT) 


ENT. 11-20 F10 XI (1) Longitudinal center of 
gravity, feet forward of 
transom. 

21-30 F10 YEL) Transverse center of gravity, 


feet to starboard. 


31-40 F10 ZI (I) Vertical center of gravity, 
feet above baseline (keel). 
fast 1-10 F10 Discrete weights termination 
card tag. 


(must be less tham 0.0) 


The maximum runber of discrete weights is 201. They are 
assumed in symmetric distribution about the longitudinal 
center plane, so that only weights on starboard side should 


be specified, with atransverse location greater than zero. 


CK 3 Craft Geometric Input 


— a s eee ee ee ee ee سیت سے سے‎ ee ee ee ee ee < ee ee — oe eee oe me pe 


No options 


Card Column Format Name Entry 
1 E’ 1 00301 
5-15 ЕТО МТА (1) Number of stations оп port 
sidewall. 
mo 25 F10 ШЕЛІГІ 22) Nunber Om stations on 
Starboard sidewall. 
26-35 F10 NSTA (3) Number of stations on bow 
seal. 


36-15 F10 NSTA (4) Number of stations on stern 
seal. 
6-5 F10 ЖЕСІ Total craft length. 
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Rote: 

The maximum number of stations is 11 for each element. The 
olan Will automatically subdivide the element into К-1 
sections (N is the number of stations) and calculate the x 


and y coordinates of each station. 


жақты»‏ = س ee‏ سس > س سے سے سے س جس سے سک سے سے 


Mption tag 1 = Sidewall input 


Card Column Format Nane Entry 
1 5 | 00401 
E15 F10 Yon y distance from centerline to 
sidewall, feet. 
16 = 25 F10 XLSW Average wetted length of 
Sidewall, feet. 
26-35 F10 CESE Leakage orfice coefficient of 
sidewall. 
36-45 PEL CDS (Os low drag coefficient ог 
sidewall. 
20-55 PTO AVBASW Average bean of  sidewall, 
feet. 
Note: 
Enger this option tag the program will also read the 


ENdewall integral table file stored on a special unit for 


mist access. 


ーー ーー ーー ーー 一 — — 


- — — سے‎ ーー ーー ーー ーー = —P w uUo uuo db GER m 
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Card Column Format Name Entry (cont) 


1 6-15 F10 XSST x coordinate of seal hinge, 


feet forward of tramsom. 


6-25 F10 ASSI 27-٦٠31٣٢٣۶ seal. binge, 
feet above keel. 

26-5 F10 ALEAK Base leakage area, sq ft. 
(ЗИ 15 nominal.) 

56-45 F10 Cees Seal leakage orifice 
مات ہ٢٦٢ ٤٣ہ وت‎ 

46-55 F10 7185 5 Length of rear support cable, 
feet. 
(1575 feet is full dorna) 

26-65 r10 DPBS Pressure differential betrean 


stern seal and plenum chamber, 
psf. 
29-75 F10 XLF Length of leading edge of the 


seal, feet. 


٣۰ © Bowseal Input 


 — s m —‏ اا = — > = Г M‏ سه چ = — د کے سے کس 


ion tag 1 = First data card 
Сага Column Format Name Entry 
1 1-5 I 00601 
6-15 F10 XBSI X Coordinate of seal Ringe, 


feet loreard of transon. 

15-25 F10 GEES Seal leakage Оисе 
coefficient. 

26-35 F10 DES Pressure differential between 
bow seal and plenum chamber, 
DSL. 

36-45 F10 1و‎ 2 coordinate of seal hinge, 


feet above keel. 
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Card Column Format Name Entry (cont) 


1 io つつ F10 ELHAXB Length of rear support cable, 
feet. 
(1.875 feet is full down.) 
565 RIO XBF Length of seal leading eäge, 
feet. 
пот 75 Е10 BLEAK Base leakage area, sq ft. 


Merron tag 2 = Second data card 


Card Column Format Name Entry 
1 1-5 I 00602 
E15 F10 CENCAB Length of middle Support 


cable, feet. 
ШЕР С-Есес is full down.) 


Option tag 1 = Plenum geometry 


Card Colunn Format Name Entry 
1 1-5 1 00701 
ст 15 F 10 LBV Plenum length at water 
surface, Ieet. 
Що - 25 F10 XBBW Plenum width at water surface, 
feet. 
5-35 F10 XPWV X coordinate of pressure wave 


pivot point, feet forward of 


+0 7. 
36-45 F10 WIDTH Pienum width at deck, feet. 
46-55 F10 XL Plenum length at deck, reet. 
267-65 F10 XCPO X coordinate of center of 


pressure, feet forward of 


transon. 
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Card Column Format Name Entry (cont) 


1 26275 F10 BUBHGT FI enum average height, feet. 


Option tag 2 = Critical Froude number 


Card Column Format Nane Entry 
1 = > 1 00702 
ре 15 FIG FNCRIT Froude number corresponding to 


hump speed. 


s, — a‏ < کے m‏ ید < w w << w s‏ > ہ ee ee‏ س جد سے ج 


No options - 


Card Column Format Name Entry 
1 09 I 00801 
б= 15 Е10 ХРО X Coordinate of propeller 
center, feet forward Qf 
5781 1501 
0 25 F10 Y pO y distance from centerline to 


propeller center, feet. 
fS 35 F10 Z PO 27 Coordinate of propeller, 


feet above Keel. 


Card Column Format Name Entry 
1 1-5 I 00901 
6-15 F 10 XRO x Or late of centroid of 
I udder, feet forward DE 
transom. 
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Card Column Format Name 
1 16-25 F10 YR 
(Dee うつ F10 ZRO 

36-45 Е10 RSPAN 

46-55 ۶10 RASPR 

56-65 Ғ10 RAREA 
66-75 F10 RIC 


Card Column Format Nane 
1 = 5 1 
0-15 FLO XLAERO 
16-25 F10 ВЕАМ 
Note: 
The reference lencth and 


nondimmensionalize the wind 


Width 


— — ーー 


y distance from centerline to 
puddcpeceoentrord, feet. 
ZEEcoomdgrnate of centroid of 
feet above keel. 


feet. 


rudder, 
Rudder span, 
Rudder aspect ratio. 
puddcrcaredq. sq ft. 


Average thickness ratio of 


rudder Section. 


Entry 
01001 
Reference length, reet. 
Reference width, feet. 


are those used to 


tunnel data used in the program. 


Block 11 Sea State Input 
Option tag 1 - Amplitude versus frecuency 
Card Column format Name Entry 
1 ШЕ с І 01101 
6-15 F10 NWAVE Number of wave components, 
maximum of 10. 
16-25 F10 BETAD Initial wave heading, degress. 
(180 degrees = ahead waves) 


32 














Mara Coluan Format 


a ғы» aÀ—À — — کے‎ 


Nane 


F10 
F10 


[Е То 
ШІ- 20 


Из, OMEGA (1) 


A9 (I) 


Option tag 2 


1 1-5 T 
6-15 F10 NwAVE 
16- 25 F10 BETAD 
EN 1-10 F10 WAVLEN (T) 
11-20 F10 AW (I) 
Note: 
En this 


of regular wave components. 


component is input. If 


ーー 
-- 


Option tag 3 


period linits 


Card Column Format Nane 
1 1-5 it 

6-15 F10 NWAVE 

16-25 F10 BETAD 

26-35 F10 SHTRV 

36-45 F10 PERL 


For a reguiar Sea case 
the 


(calm water), only the first card 


Entry (cont) 


Wave frequency, rad/sec. 


wave amplitude, feet. 


Amplitude versus wave length 


Entry 
01102 


of 


maximum of 10. 


Number wave components, 


Initial wave headina, degrees. 
(180 degrees - ahead waves) 
Wave length, feet. 


Wave amplitude, feet. 


options irregular seas are constructed by addition 


only one 
number of components is Zero 


may appear. 


Fave Components by calculation for given 


Entry 
01105 


Gf wave 


max1mum of 10. 


Number components, 


Initial wave heading, degrees. 


ahead waves) 
ОГ ВУЗ +٤ 


(180 degrees 
height 
Teet. 


Shortest significant period of 


AVI ge 


vaves, 


waves, seconds. 
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ane Entry (cont) 
1 2-55 F 10 PER Longest Significant period gf 


waves, seconds. 


Option tag 4 = Wave components by calcuiation for given 


frequency limits 


Card Column Format Name Entry 
1 1-5 m 01104 
5-15 F10 NWAVE Number of wave components, 


Ша КОШО лл от 10. 


5-25 F10 BETAD Initial vave heading, degrees. 
(180 degrees - ahead waves) 
5-35 Ғ10 SHTWV Average height of 1/3 highast 


waves, feet. 

36-45 FIO WNN Lowest significant frequency 
OÍ waves, Hertz. 

46-55 F10 WWX Highest significant frequency 


of waves, Hertz. 


= 


In options 3 and 4 the wave components are calculated in 
Order to ease the input requirements for irregular sea 


Cases. 


Beck 12 Initial Conditions 


سے ーー ーー ーー ーー ーー ーー‏ سے سی 


Card Column format Name Entry 
1 1-5 T 01201 
6715 F10 UO initial speed, knots. 
16-25 F10 THETO NUNC pitch angie (positive 


Senese sp, in degrees). 
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Card Column Format Name Entry (cont) 


1 шъ 35 PTO DSO 7۰۰٣۰۷۰٦٦٦47۹3 ٣٢٣ج‎ center of 
gravity, inches. 
36-45 F10 DELPI initial bubble gage pressure, 
psf. 


Block 13 End of Case 


a a ہہ سس :+ سے سس ست سس جک سے سے سے س‎ шы Do سسس‎ ee 


Card Column Format Nane 77 
1 1-5 I 013 
Note: 


This card sianals the end of data for a given case. Two 
consecutive 013 cards with no data between vili stop the job 


unduly. 


Block 14 rnd oz Run 


No options 


Card Column Format Name Entry 
1 1-5 1] 014 
Note: 


This card signals the end of the run. this will be indicated 
by the message: COMPLETED ALL RUNS. 


ーー ーー ーー ーー ビー ーー ーー ee ےک‎ ne 


Option tag 1 = Starboard thrust input 
ormat Name Entry 


1 (025 1 01601 
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Card Coiumn 
1 6-15 
16-25 

26-35 

2,3,. 1-80 
1-80 


mote: NPS must not 


Option tag 2 


6-15 
Шо 25 
0-35 
ENS. 1-80 
1-80 


EIO 
ВЕТО 


= Port 


Fornat 
1 


F 10 
F10 
F10 
ВЕТО 
ВЕТО 


Name 


THST1 


NPS 


STHS 
TIS (J) 


THSTS (3) 


HAU “(eon } 


Value of the thrust for the 
Starboard screw. This is an 
estimate or kept fixed 
depending on option switch in 
BEOCK. 00105. 

Number ог data points in the 
Пи nap: Left plank iftwano 
thrust map is used. 

Sidefthr ust coefficient. 

Data points for independent 
variable time in ascending 
order. 

Corresponding data points for 


the dependent variable. 


Mental ge nay25. 


йе е еше 


THST2 
NPP 
STHP 
TIP (J) 
THSTP (J) 


Entry 
01602 


A RADOVE TOY Port Seren. 
ASADO port screw. 
NS зпоме tor port. screWw. 
Дещароте tor port Screw. 


As above for port screw. 


Note:  NPP must not be larger than 25. 


Option tag 3 


- Rudder motion input 
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Card Column Format Name ЕПБЕШИ (с сит) 


1 ше 15 F10 DELR Value of rudder if constant 
rudder is desired (in 
degrees). 

5-25 F10 NPR Number of data points of the 
rudder ma. Left blank 1f 
constant rudder is desired. 

ЕЕ). 1-80 ط8‎ ۱۷۵ TIR (J) Values of independent variable 
time in ascending order. 

00 8F10 DELRUD(J) Corresponding values of rudder 
angle (in degrees). 

Note: Positive rudder angle 


i TIGR rudder. 


NPR must not be larger than 25. 


Option tag 1 = Bow seal 
Card Column Format Name Entry 
1 5 T. 01701 
i15 F10 NB МОО ت۰۱‎ 7 data 001018 6 


bow seal. 


E. 1-30 OEC EE pnr) Tine points in ascending 
order. | 
EU. 1-80 ООО) Corresponding differences in 


bow seal pressures. 


Option tag 2 = Stern seal 
Card Column Format Mane entry 


1 0-5 it 01702 
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Card Column Format Name 

1 E15 F10 NS 
EN. 1-80 SINUS ІМ 
7. 1-80 ӨНІП DETS (T) 
Note: 
Ше constant pressures 


one seal alone 


we ic Title Card 


os < Á Á ee ee ымы» 


Block 1° Fan M 


سوہ — ee ee ee  —‏ مت سے سے 


Option tag 1 = 


TOO. 


aps 


BOW 


Card Column Format 


و09 1 


8-15 
16-25 
26-5 


T 


F10 
F10 
F10 


EBD Ut in block 5 and/or block 


paz 
12) 
13 
ID 


mul re 


seal fans 


Nane 


ENBFAN 
BRPM 
NPTSB 


are 


Entry (cont) 


Number of data points for the 
stern seal. 

ALO Option 1, for stern 
seal. 

As for option 1, for the stern 


seal. 


desired, those values shall] be 


DNhnJgceblock^can be used For 


Entry 
01801 


Alphanumneric data tO be 
printed as title in input 


summary. 


Kumper of bow seal fans. 
Speed of bow seal fans, rpm. 
Number of data points. 


(naximum of 25) 
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— — aa | یو‎ —— ーー ーー ーー ーー ーー C— mM — ч че» 


2,3,. 1-80 8F10 


EU. 1-80 8F 10 


Name 


READIN 


PBFAN (3) 


QBFAN (J) 


Option tag 2 = Cushion fans 


Card Column Fornat 
1 1-5 I 

6-15 F10 
16-25 F10 
26-35 E19 
36-45 F10 

EL3,. 1-80 8F10 

6,7,. 1-80 8F10 


Option tag 3 = Stern seal fans 


1 Ш- 5 І 
5-15 F10 
5:25 TUO 


Name 


ENHFAN 
EAS PBI 
۳200 


READIN 


PUFEN (J) 


QMFAN (J) 


Nane 


БМ5ҒАМ 
SKPM 


Entry (cont) 


sis a 1.0 then the 
following data points are 
read. If left blank no further 
read takes place. 

Tabular values of fan pressure 
бз erential, psf. 

Corresponding volumetric flow 


БН ЕЕЕ сес. 


aber zof’eushıon fans. 

БЕЛ ШОТ cushion fans, rpm. 
Number of data points. 
(maximum of 25) 

Мәтел а 1:20. then tae 
following data points are 
read: f lert blank no further 
read takes place. 
DapmulasSvwarues or fan pressure 
differential, psf. 
Corresponding volumetric flow 


т со ftt/sec. 


Number of stern seal fans. 


Speed of stern seal fans, rpm. 
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Card Column Format Name Entry (cont) 
1 0-5 F 10 NPTSS Number of data points. 
(naximum of 25) 
96-15 F10 READIN ОО И << а 1.0 then the 
following data points are 


۶۰۰٠۰۰۰۰٢ tere blank no further 
read takes place. 

ENS. 1-80 ВЕТО PSFAN(J) Tabular values of fan pressure 
@igeerential, psf. 

E 1-80 ВЕЛО QOSFAN(J) Corresponding volumetric flow 


Pace yen 1 t/ Sec. 


The tabular pressure data for ali the fans must be stored in 


EES Ot increasing pressure. The data is for a single fan. 





Card Column Format Name Ente) 
1 1-3 020 
2 = 0 1011 NCURV үе "Пер ой this curve on 


tnescorressponding azaph, 1.6. 
column 1 pertains to graph 1, 
БӘЛІШ 2 Со Graph 2, ге. тіс: 

0 


E S т< от more than one 


only one curve per grapa. 


curves per graph. 
2 = intermediate of three 
curves per graph. 
DE erst Of more than опе 


curves per graph. 
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Bote: Haximum number o£ curves per graph = 3. 


Maximum number of graphs = 10. 


Block 21 Tabular Summaries 


due ニー ニー ーー i i i eee ee ee 


Card Column Format Name Entry 
1 [ЕЗ ІЗ 021 
2 = 16 812 1 71 Numbers corresponding to the 


variable names to be printed 
out in summary one. 

3 1-16 812 ISJ“? Number corresponding to the 
variable names to be printeä 
out in summary two. Summary 
two ls printed out in 
successive runs without the 


cards being Tepeated. 


u Q s eee ce A ee ee aman ет vie cA MO جم‎ A A سے _— سے‎ 


Card Column Format Name Entry 
1 03 та 022 
2 ШЕ 2 I2 NGRAF Number ОЕ graphs desired, 


٣۰۰٠۰٠٠٢٢٠٢ ں٦٣:‎ ۶1 م39‎ up to and 


including 10. 


3-4 ша NDRW 00 - print plot 
01 = CALCOMP plot 
9 080۸02 26 RAYS Numbers corresponding to the 


variable names desired fcr tne 
X and Y axes of the graphs in 


alternating manner. 
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Card Column Format Nane Entry (cont) 


4 1-48 бАВ TICRD Alpha-numerical title on the 
Far St Veneer OL the CALCOHP 
٠۰۰۰۰٠٠٦۴ f print plot 15 


desired. 


— q —— —  — — — D 


In order to maintain the pitch angle curve of figure E 
the present computer model has a speed dependent shitt 
routine of XCPO of the following form: 

XCPU=XCP +A* (U*0.5921-30.0) **2+B 
After progran modifications the coefficients Á and B may be 
Mecaliculated by setting ХСРИ eguai to XCP and using 
aifferent input values for YCPO at 10 knots and 30 knots 
until the desired pitch angles are obtained. Figure 8 gives 
sone initial estimates. From the two data points the values 
for A and B can be determined replacing XCPO by its 6+07 


value again. 
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مس ape aap‏ — — ہے 








ー ーーーーーー 


020 025 





SCALE=5.00E+00 UNITS INCH. 
SCALE=S.00€-01 UNITS INCH. | 
BUGRAM 2, 1. STRAIGHT RUN AT 20 KNOTS 


-OT IS PITCH ANGLE VERSUS TIME 


103 













LE=5.00E+00 UNITS INCH. 104 
ILE=1. QOE+00 ШИЛІ БЫ. ғ. 
RAM 2, 1. STRAIGHT RUN AT 20 KNOTS 
| M JISPLBUEMEN I UERSUS TIME 












5.00E«00 UNITS INCH. 

E=5.00€+00 UNITS INCH. | 

IHM 2, 1. STRAIGHT RUN AT 20 KNOTS 
15 PLENUM PRESSURE VERSUS TIME 
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РАО А 


Lr а ーーmーーーーー ON A ee 








— 0. 








A A —— — 


_ 9685 - 010 mi ci — Pd ae 025 020 03 





хо. 00Е+00 UNITS INCH. 

5. 00Е+01 UNITS INCH. 

AM 2, 1. STRAIGHT RUN RT 20 KNOTS 
IS THRUST STARBOARD VERSUS TIME 
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 ——— 


010 25 (25 925 «ЖЕГІ 






5. D0E+00 UNITS INCH. 
F=2.006-0! UNITS INCH. 


AM 2, 2. STRAIGHT RUN AT 20 KNOTS 
IS PITCH ANGLE VERSUS TIME 


] 
RE 
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ーー = = = 
AAA — ا ص‎ on 


A 31 02: 







.=3. 00E+00 UNITS INCH. 108 
=1.00E+00 UNITS INCH. | E 
aM 2, 2. STRAIGHT RUN AT 20 KNOTS 
Bc DISPLACEMENT ШЫ ٦ 








— = — م یس مہ‎ — АН 
(ao 
rm 


019 9:0 дес 025 930 





:5.00E«00 UNITS INCH. 
#5. 00E+00 UNITS INCH. 


AM 2, 2. STRAIGHT RUN AT 20 KNOTS 
[S PLENUM PRESSURE VERSUS TIME 


109 








5. 008+00 UNITS INCH. 
L-o.QOE«01 UNITS INCH. 


AM 2, 2. STRAIGHT RUN AT 20 KNOTS 
1S THRUST STARBOARD VERSUS TIME 


110 











010 015 Ё 77 920 | O25 3080 = р 035 





Е-5.00Е+00 UNITS INCH. 
=5,00E-01 UNITS INCH. 


AM 3. 1. STRAIGHT RUN AT 20 KNOTS 
IS PITCH ANGLE VERSUS TIME 


ILI 





ОСАГЕ 5, 00Е+00 UNITS INCH. 


SCALE=1.00E+00 UNITS INCH. 
AUGRAM 3, 1. STRAIGHT RUN AT 20 KNOIS 
EE 15 7 DISPLACEMENT UERSUS TIME 





RLSP IL 











9. 00E+00 UNITS INCH. 
то. О00Е+00 UNITS INCH. 


aH 3, 1: STRAIGHT RUN AT 20 KNOTS 
IS PLENUM PRESSURE VERSUS TIME 


TES 





| PEOT 12 


ーーーーーーー - 






«5.00E«00 UNITS INCH. 
S5. 00E«01 UNITS INCH. 


AM 3. 1. STRAIGHT RUN AT 20 KNOTS 
IS THRUST STRRBORRD UERSUS TIME 


PEA 





TT 15 


gH `; 





010 035 022 229 030 Ua 






=5.00E+00 UNITS INCH. 
E=1.00E-01 UNITS INCH. 


RAM 3, 2. STRAIGHT RUN AT 20 KNOTS 
| 15 BITCH ANGLE UERSUS TIME 


ТБ 





PLOT 14 





010 245 020 025 039 057 





о UNITS INCH. | 
.QOE«00 UNITS INCH. 


۴ M 3. 2. STRAIGHT RUN AT 20 KNOTS 
IS 2 DISPLACEMENT VERSUS TIME 


116 | 





029 = 225 7703 






5.00E«00 UNITS INCH. 

E 00€+00 UNITS INCH. 
AM 3- 2. STRAIGHT RUN AT 20 ص50‎ 
IS PLENUM PRESSURE VERSUS TIME 


117 


| 





PLOT 16 





a 
- 


0:0 05 02Q 075 039 уж. 






5.00E+00 UNITS INCH. 
=5.00E+01 UNITS INCH. | 
ІНМ 3, 2. STRAIGHT RUN RT 20 KNOTS 

| 15 THRUST STARBOARD VERSUS TIME 
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PLOT 17 











E=5.00€+00 UNITS INCH. 
LE=5. 00-01 UNITS INCH. 


RAM 2. TURN AT 20 KNOTS WITHOUT A.D. 
.]S PITCH ANGLE UERSUS TIME 
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DEC alee) 








Е«Б.00Е-00 UNITS INCH. 
Ez2.00F+00 UNITS INCH. 


ARM 2. TURN AT 20 KNOTS WITHOUT. 
IS PITCH RATE VERSUS TIME 


120 


H 。 


p 

















PLOT 19 








010 72 023 025 Е, 05 






0, 00E+00 UNITS INCH. 

2. 00Е+00 BEES us i 

AM 2, TURN RT 20 KNOTS HITHOUT R.D. 
15 ROLL ANGLE VERSUS TIME 
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PLOT 20 





=5.09E+009 UNITS INCH. 

6. 00Е+00 UNITS INCH. 

2, TURN AT 20 KNOTS WITHOUT R.D. 
ROLL RATE VERSUS TIME 


RAM 
BIS 


0 











5 
[ 


:5.00E«00 UNITS INCH. 
15. 00Е-01 UNITS INCH. 


AM 3, TURN AT 20 KNOTS WITHOUT RD: 
IS PITCH ANGLE VERSUS TIME 


123 


010 015 020 025 o 039 O < 





LOT ص2‎ 





SCALE=5.00E+00 UNITS INCH. 
SCALE=1.00E+00 UNITS INCH. ーー _- 
JÓBRRM 3. TURN AT 20 KNOTS WITHOUT 1۱ 
Т 15 PITCH RATE VERSUS TIME 
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GARE =5.00£+00 UNIT TS INCH. 
GALE=2.008+00 UNITS INCH. 


GRAM 3. TURN AT 
JT IS ROLL ANGLE 


i - 


031 


[ 3 


20 KNOTS WITHOUT Re 


VERSUS TIME 


125 





LE=5. 00E+00 UNITS INCH. 
-5.00E+00 UNITS INCH. 


AM 3, TURN AT 20 KNOTS WITHOUT R.D- 
IS — ROLL RATE UERSUS TIME 


126 





EM LLL LL 


PLOT 25 











010 gu 022 225 oo بے‎ 02“ 





F«5.00£«00 UNITS INCH. 

-5. 006-01 UNITS INCH. | 
ШОУ, TURN AT 20 KNOTS With R.D. 
LO PITCH ANGLE VERSUS TIME 


Po 


121 








CALE=S.00F+00 UNITS INCH. 

GALE=1.006+00 UNITS INCH. 

۰۰٦ 2, TURN AT 20 KNOTS WITH R.D. 1 
9 9 PITCH RATE VERSUS TIME 


128 





PLOT 27 





92% 925 930 


Ш ИС -5.00ғ-00 UNITS INCH. 
SOALE=1. 00E+00 UNITS INCH. | | 
¡OGRAM 2, TURN AT 20 KNOTS WITH A.D. 1 
DI IS ROLL RNGLE VERSUS TIN 


o) 
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m af? a = سے سی ہے‎ 








SCALF-5.00E+00 UNITS INCH. 
30RLE--2.00£«00 UNITS INCH. 
WGRAM 2, TURN AT 20 KNOTS ٢ ۰ا‎ 
Ж Іс ROLL RATE VERSUS ТІЗЕ 


130 





CRLE-5.00E«00 UNITS INCH. 
CALE=2.00E-01 UNITS INCH. 
OGRAM 3. 


ОГ IS PITCH ANGLE 


131 


енер олар ممت میت‎ s 


— — — 
> 


— ےج ہے ےچ ا ححے مہ ےم‎  —n 
= 


TURN AT 20 KNOTS WITH RU, 


VERSUS TIME 





| 





PLOT 30 








SCALE=5. 00E+00 UNITS INCH. 

SCALE-=1.00E+00 UNITS INCH. 

WERAM 3. TURN AT 20 KNOTS WITH R.D. 1 
ШІ І5 = PITCH RATE VERSUS TIME 





PLOT 31 





—- — o ... 


UCS 010 915 020 ن02‎ 030 


SWHLE=S。0UE+00 UNITS INCH. 
SCHLE=1.00E+00 UNITS INCH, 


106081 3 TURN AT 20 KNOTS WITH R.D. 1 
ОТ 1S ROLL ANGLE VERSUS TIME 


155 





PED 32 


し レー ミー O A ee s.s 





SCALE=5.00E+00 UNITS INCH. 

SCALE=5.00E+00 UNITS INCH. 
JØGRAM 3. TURN AT 20 KNOTS WITH R.D. 1 
ОТ 15 = ROLL RATE VERSUS TIME 











PLOT 33 


= ee 


—-—r a 
— F" 
—  —— U... U U U. 
— 


A en s sa 
а 
G: 


SCALE=5.00E+00 UNITS INCH. 
SCALE=2.00E-01 UNITS INCH. 

ШОВАм 3, TURN AT 20 KNOTS WITH R.D. 
Me 1S PITCH ANGLE VERSUS TIME 


| 
1 


| 


195 





PLOT 34 


єз 
со 
єз 





dl nn 


-E=9.00E+00 UNITS INCH. 
ша. ПОЕ+00 UNITS INCH. 





کہ 


2 


VERSUS Tino 


TURN AT 20 KNOTS WITH RD. 
РІТОН RATE 


) H Е M = と 
E 


136 





we ہے‎ > A 





XRLE-5.00E«00 UNITS INCH. 
SCALE=1.00E+00 UNITS INCH. 


¡GRAM 3. 


TURN AT 20 KNOTS WITH RD: 


OT 15 ROLL ANGLE VERSUS TIME 


d 





PLOT 36 


mga, Àm sb سے‎ m q = 








SCALE=5.00E+00 UNITS INCH. 
SCHLE=5.00E+00 UNITS INCH. 


ROGRAM 3. TURN AT 20 KNOTS HITH R.D. Z 
LOT IS X ROLL RATE VERSUS TIME 


| 
| 
138 
| 








SCALE-5.00E+00 UNITS INCH. 
SCALE=2.00E-01 UNITS INCH. 
IJBHRM 3. TURN 8۲٢ 20 KNOTS HITH R.D 
Mm IS PITCH ANGLE VERSUS IIML 


139 





PLOT 38 





EN —5.00E-00 UNITS INCH. 
ши 1.00E-00 UNITS INCH. 


JGRAM З, TURN AT 20 KNOTS WITH 8.D. 3 
Mm 1S PITCH RATE VERSUS TIME 





PLOT 39 


| VA 


Алл 


へ 、 











CALE=5.006+00 UNITS INCH. 

CALE=1. 00€+09 UNITS INCH. 

JGRAM 3. TURN AT 20 KNOTS WITH R.D. 3 
JT IS ROLL ANGLE | VERSUS TIME 





| 
I 





E=5.00E+00 UNITS INCH. 
E=5.00E+00 UNITS INCH. 


RAM 3. TURN AT 20 KNOTS WITH P.D. 3 
| IS ROLL RATE VERSUS TIME 
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CALE:5.00£+00 UNITS INCH. 
CALE =2. 006-01 UNITS INCH. 

UJGRAM 3. TURN AT 20 KNOTS WITH R. 
JI 1S PITCH ANGLE VERSUS ІП 


d 
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PLOT 42 





۷۰۰۰۰۰7 ٥0+00 UNITS INCH. 
Cart =1.00E+90 UNITS INCH. 


ӘбААМ 3. TURN BT 20 KNOTS WITH R.D. 4 
Jl IS PITCH RATE pU» ۷٣۳ 
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GALE=S.00£+00 UNITS INCH. 

OBLE=:1. 00£+00 UNITS INCH. 

UGRAM 3, TURN ПТ 20 KNOTS WITH &.0. 4 
JI IS ROLL ANGLE VERSUS TIME 
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PLOT 44 











JUBLE-5.Q0E«00 UNITS INCH. 

BEALE=2.00E+09 UNITS INCH, | | 
ШӘНГН 3. ТОҢЫ 81 20 KNOTS HITH R.D. d 
m 15 ROLL RATE VERSUS ГІНЕ 
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٦۰5 00۲+00 UNITS INCH. 
25 —2.00E-01 UNITS INCH, 

٦ / | г 
¡UGRAM 3, TURN ET 20 
OT 1S PITCH ANGLE 
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EESE-5.00F«00 UNITS INCH: 

ШИС-1.00:-00 UNITS INCH. | š 
GRAM 3. TURN RT 20 KNOTS HITH R.D. 5 
ОТ 15 PITCH RATE VERSUS TIME 
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9:9 2-> Des 020 9sQ 93 









F«5.00£400 UNITS INCH. 
۰۱۰000+00 UNITS INCH. 


ЗАМ 3. TURN 91 20 KNOTS HITH 8.0. 5 
1S ROLL ANGLE VERSUS TIME 
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ЕНЕ шш > 


ÍA 
ГА V 


552 


8 : A D o oM е m 
, | 7 nd 253 Ce vd 929 Em 
J 





ШЕ + 00 UNITS INCH: 

5. п0Е+02 UNITS INCH. ーー 
AN 3. TURN AT 20 KNOTS HITH B.D. + 
IS AGLL RATE VERSUS TIME 








| ーーーーーーーー 一 ーーー my a 
ہچ‎ У ` г с,” Use + б 1. 
| [VIVIS] ve J “еч vo Шо 


8۰۰۲۰۰00۷۲۰00 UNITS INCH. 

NEU. 00E-01 UNITS INCH. 

0+001 3: TURN AT 20 KNOTS ИТІН ч.0. b 
Ul IS PITCH ANGLE VERSUS TIME 


ku 
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PLOT 50 
| 


5668 3, TURN AT 20 KNOTS HITH R.D. 
UT IS PITCH RATE UERSUS TIME 


C 


۱ 
XBEE-5.00E«-00 UNITS INCH. 
IBHLE-1.Q00E«00 UNITS INCH. 











BLE-5.QOE+09 UNITS INCH. 
JOBEE-:1.00E«00 UNITS INCH. 


IJGRAM 3, TURN AT 
BI 1S ROLL ANGLE 


av 
DE 
20 KNOTS WITH R.D. 


— سے‎ ーー 


VERSUS TIME 
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5. 00E+00 UNITS INCH. 

5. 00E+00 UNITS INCH. 

RAM 3. TURN AT 20 KNOTS WITH B.D. € 
Іс ROLL RATE VERSUS TIME 











DE —5.00E«O00 UNITS INCH. 

CALE=2.00E-01 UNITS INCH. | | 

ØGRAM 3. TURN AT 20 KNOTS WITH A.D. 7 
OT 1S PITCH ANGLE VERSUS TIME 
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SCALE=5.00E+00 UNITS INCH. 
SCALE=1:.00E+00 UNITS INCH. 


WGRAM 3, TURN AT 20 KNOTS: WITH R-N. 
ща 15 PITCH RATE VERSUS TIME 
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| 
| 25 010 0:5 023 025 — 030 


| 
| | 


۲۰۰۰٦-١٠0٢٣٣ +00 UNITS INCH. 
pSHEE-1.00E«00 UNITS INCH. 
OGROM 3- TURN AT 20 KNOTS WITH P.D. 
(OT 1S RØLL ANGLE VERSUS TIME 
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PLOT 56 








3CRLE-5. 006.00 UNITS INCH. 

PORLE-5.00€«00 UNITS INCH. 

UGRAM 3, TURN AT 20 KNOTS WITH R.D. 7 
E |^ RULL RATE VERSUS TIME 





PLOT 57 


JURLEs5.00E«00 UNITS INCH. 

SCALE=2.00E-01 UNITS INCH. | 
WCRAM 3. TURN AT 20 KNOTS HITH R.D. G 
= 15 PITCH ANGLE VERSUS TIME 








! | | | A 1 gg 7 


SCALE=5.00E+00 UNITS INCH. 
PURALE==1.00€+00 UNITS INCH. 


IUBBHM 3. TURN AT 20 KNOTS WITH R.D. € 
ШІ5 PITCH RATE VERSUS TIME 
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my to s _ a J ーー 





) -:5.00E+00 UNITS Nun: 161 

UNITS INCH.‏ 7 >> تت 

GRAM 3. TURN AT 20 KNOTS WITH ۰ 8 
OT IS ROLL ANGLE VERSUS TIME 


ге! 








CALE=5.00E+00 UNITS INCH. 

SCALE=S.00£+00 UNITS INCH. | 

GRAM 3. TURN AT 20 KNOTS WITH R.D. 6 
E I5 ROLL RATE VERSUS TIMt 
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